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We provide a brief review of recent progress in the study of mesons using QCD's Dyson- 
Schwinger equations. Along the way we touch on aspects of confinement and dynamical chiral 
^ • symmetry breaking but in the main focus upon: exact results for pseudoscalar mesons, in- 

^ I eluding aspects of the 77-77' problem; a realisation that the so-called vacuum condensates 

n ^ . are actually an intrinsic, localised property of hadrons; an essentially nonperturbative pro- 

cedure for constructing a symmetry-preserving Bethe-Salpeter kernel, which has enabled a 
0\J ■ demonstration that dressed-quarks possess momentum-dependent anomalous chromo- and 

^ I electromagnetic moments that are large at infrared momenta, and resolution of a longstand- 

^L/ ' ing problem in understanding the mass-splitting between p- and ai-mesons such that they 

^— V I are now readily seen to be parity partners in the meson spectrum; features of electromagnetic 

^^ ■ form factors connected with charged and neutral pions; and computation and explanation of 

I ■ I valence-quark distribution functions in pseudoscalar mesons. We argue that in solving QCD, 

f— ^ ■ a constructive feedback between theory and extant and forthcoming experiments will enable 

constraints to be placed on the infrared behaviour of QCD's /3-function, the nonperturbative 
quantity at the core of hadron physics. 
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I. INTRODUCTION 

A hundred years and more of fundamental research in atomic and nuclear physics 
has shown that ah matter is corpuscular, with the atoms that comprise us, themselves 
containing a dense nuclear core. This core is composed of protons and neutrons, referred to 
collectively as nucleons, which are members of a broader class of femtometre-scale particles, 
called hadrons. In working toward an understanding of hadrons, we have discovered that 
they are complicated bound-states of quarks and gluons. These quarks and gluons are 
elementary, pointlike excitations, whose interactions are described by a Poincare invariant 
quantum non-Abelian gauge field theory; namely, quantum chromodynamics (QCD). The 
goal of hadron physics is the provision of a quantitative explanation of the properties of 
hadrons through a solution of QCD. 
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Quantum chromodynamics is the strong-interaction part of the Standard Model of 
Particle Physics and solving QCD presents a fundamental problem that is unique in the 
history of science. Never before have we been confronted by a theory whose elementary 
excitations are not those degrees-of-freedom readily accessible via experiment; i.e., whose 
elementary excitations are confined. Moreover, there are numerous reasons to believe that 
QCD generates forces which are so strong that less-than 2% of a nucleon's mass can be 
attributed to the so-called current-quark masses that appear in QCD's Lagrangian; viz., 
forces capable of generating mass from nothing, a phenomenon known as dynamical chiral 
symmetry breaking (DCSB). 

Neither confinement nor DCSB is apparent in QCD's Lagrangian and yet they play 
the dominant role in determining the observable characteristics of real-world QCD. The 
physics of hadrons is ruled by emergent phenomena such as these, which can only be 
elucidated through the use of nonperturbative methods in quantum field theory. This is 
both the greatest novelty and the greatest challenge within the Standard Model. We must 
find essentially new ways and means to explain precisely via mathematics the observable 
content of QCD. 

The complex of Dyson-Schwinger equations (DSEs) is a powerful tool, which has 
been employed with marked success to study confinement and DCSB, and their impact on 
hadron observables. This will be emphasised and exemplified in this concise review, which 
describes selected recent prog ress in the study of mesons. It complements and extends 
earlier and other efforts [l|-[lO|. 



II. HADRON PHYSICS 

The basic problem of hadron physics is to solve QCD. This inspiring goal will only be 
achieved through a joint effort from experiment and theory because it is the feedback be- 
tween them that leads most rapidly to improvements in understanding. The hadron physics 
community now has a range of major facilities that are accumulating data, of unprecedented 
accuracy and precision, which pose important challenges for theory. The opportunities for 
researchers in hadron physics promise to expand with the use of extant accelerators, and 
upgraded and new machines and detectors that will appear on a five-to-ten-year time-scale, 
in China, Germany, Japan, Switzerland and the USA. A short list of facilities may readily 
be compiled: Beijing's electron-positron collider; in Germany - COSY (Jiilich Cooler Syn- 
chrotron), ELS A (Bonn Electron Stretcher and Accelerator), MAMI (Mainz Microtron), 
and FAIR (Facility for Antiproton and Ion Research) under construction near Darmstadt; 
in Japan - J-PARC (Japan Proton Accelerator Research Complex) under construction 
in Tokai-Mura, 150km NE of Tokyo, and KEK, Tsukuba; in Switzerland, the ALICE and 
COMPASS detectors at CERN; and in the USA, both the Thomas Jefferson National Accel- 
erator Facility (JLab), currently being upgraded, with new generation experiments expected 
in 2016, and RHIC (Relativistic Heavy Ion Collider) at Brookhaven National Laboratory. 

Asymptotic coloured states have not been observed, but is it a cardinal fact that 
they cannot? No solution to QCD will be complete if it does not explain confinement. 
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This means confinement in the real world, which contains quarks with light current-quark 
masses. This is distinct from the artificial universe of pure- gauge QCD without dynamical 
quarks, studies of which tend merely to focus on achieving an area law for a Wilson loop 
and hence are irrelevant to the question of light-quark confinement. 

In stepping toward an answer to the question of confinement, it will likely be necessary 
to map out the long-range behaviour of the interaction between light-quarks; namely, QCD's 
^-function at infrared momenta. In this connection it is noteworthy that the spectrum of 
meson and baryon excited states, and hadron elastic and transition form factors provide 
unique information about the long-range interaction between light-quarks and, in addition, 
the distribution of a hadron's characterising properties - such as mass and momentum, 
linear and angular - amongst its QCD constituents. The upgraded and promised future 
facilities will provide data that should guide the charting process. However, to make full 
use of that data, it will be necessary to have Poincare covariant theoretical tools that enable 
the reliable study of hadrons in the mass range 1-2 GeV. Crucially, on this domain both 
confinement and DCSB are germane. 

It is known that DCSB; namely, the generation of mass from nothing, does take place 
in QCD. It arises primarily because a dense cloud of gluons comes to clothe a low-momentum 



quark 1]|. This is readily seen by solving the DSE for the dressed-quark propagator; i.e., 
the gap equation, which yields the result illustrated in Fig.[TJ (In our Euclidean metric: 
{1/1,1^} = 26^^,; -il = 7/,; 75 = 747i7273, tr [747^7^,7^7^] = -Ae^yf„-, a^u = (i/2)[7^,7,,]; 
a ■ b = y^i=i^i^i'i ^'^d -f/i timelike =^ P^ < 0. More information is available in App. A 
of Ref. [5|.) However, the origin of the interaction strength at infrared momenta, which 
guarantees DCSB through the gap equation, is currently unknown. This relationship ties 
confinement to DCSB. The reality of DCSB means the Higgs mechanism is largely irrelevant 
to the bulk of normal matter in the Universe. Instead the most important mass generating 
mechanism for light-quark hadrons is the strong interaction effect of DCSB; e.g., one may 



identify it as being responsible for 98% of a proton's mass [l5|, |16|. There is a caveat; 
namely, as so often, the pion is exceptional. Its mass is given by the simple product of 
two terms, one of which is the ratio of two order parameters for DCSB whilst the other is 
determined by the current-quark mass (Sec. lVV-ip . Hence the pion would be massless in 
the absence of a mechanism that can generate a current-mass for at least one light-quark. 
The impact of a massless, strongly-interacting particle on the physics of the Universe would 
be dramatic. 

It is natural to ask whether the connection between confinement and DCSB is acci- 
dental or causal. There are models with DCSB but not confinement, however, a model with 
confinement but lacking DCSB has not yet been identified (see, e.g., Sees. 2.1 and 2.2 of 
Ref. [tI]). This leads to a conjecture that DCSB is a necessary consequence of confinement. 
It is interesting that there are numerous models and theories which exhibit both confine- 
ment and DCSB, and possess an external control parameter such that deconfinement and 
chiral symmetry restoration occur simultaneously at some critical value of this parameter; 



e.g., quantum electrodynamics in three dimensions with Nf electrons 17Hl9l|. and models 
of QCD at nonzero temperature and chemical potential 2Ch28|]. Whether this simultane- 
ity is a property possessed by QCD, and/or some broader class of theories, in response to 
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Rapid acquisition of mass is 
^effect of gluon cloud 



— m = (Chiral limit) 

— m = 30MeV 

— m = 70 MeV 




p [GeV] 

FIG. 1: Dressed-quark mass function, M{p) in Eq. ([2|): solid curves - DSE results, explained in 
12, ll3j, "data" - numerical simulations of lattice-QCD [IJ]. (NB. m = 70 MeV is the uppermost 
curve and current-quark mass decreases from top to bottom.) One observes the current-quark 
of perturbative QCD evolving into a constituent-quark as its momentum becomes smaller. The 
constituent-quark mass arises from a cloud of low-momentum gluons attaching themselves to the 
current-quark. This is dynamical chiral symmetry breaking (DCSB): an essentially nonperturbative 
effect that generates a quark mass from nothing; namely, it occurs even in the chiral limit. (Figure 
adapted from Ref. [ll[.) 



changes in: the number of light-quark flavours; temperature; or chemical potential, is a 
longstanding question. 

The momentum-dependence of the quark mass, illustrated in Fig.[Tl is an essentially 
quantum field theoretic effect, unrealisable in quantum mechanics, and a fundamental fea- 
ture of QCD. This single curve connects the infrared and ultraviolet regimes of the theory, 
and establishes that the constituent-quark and current-quark masses are simply two con- 
nected points separated by a large momentum interval. The curve shows that QCD's 
dressed-quark behaves as a constituent-quark, a current-quark, or something in between, 
depending on the momentum of the probe which explores the bound-state containing the 
dressed-quark. It follows that calculations addressing momentum transfers Q^ > M^, 
where M is the mass of the hadron involved, require a Poincare-covariant approach that 
can veraciously realise quantum field theoretical effects [29t] . Owing to the vector-exchange 
character of QCD, covariance also guarantees the existence of nonzero quark orbital angular 
momentum in a hadron's rest-frame 



30-132]. 



The dressed-quark mass function has a remarkable capacity to correlate and to con- 
tribute significantly in explaining a wide range of diverse phenomena. This brings urgency 
to the need to understand the relationship between parton properties in the light-front 
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frame, whose peculiar properties simplify some theoretical analyses, and the structure of 
hadrons as measured in the rest frame or other smoothly related frames. This is a problem 
because, e.g., DCSB, an established keystone of low-energy QCD, has not explicitly been 
realised in the light-front formulation. The obstacle is the constraint k^ := k^ + k^ > 
for massive quanta on the light front [33] • It is therefore impossible to make zero momen- 
tum Fock states that contain particles and hence the vacuum is "trivial". On the other 
hand, it is conceivable that DCSB is inextricably tied with the formation and structure of 
Goldstone modes and not otherwise a measurable property of the vacuum. This conjec- 



ture is being explored |34l436l| and is something about which more will be written herein 
(Sec. lVV-2p . In addition, parton distribution functions, which have a probability interpre- 
tation in the infinite momentum frame, must be calculated in order to comprehend their 
content: parametrisation is insufficient. It would be very interesting to know, e.g., how^if 
at all, the distribution functions of a Goldstone mode differ from those of other hadrons 



III. CONFINEMENT 

It is worth stating plainly that the potential between infinitely-heavy quarks measured 
in numerical simulations of quenched lattice-regularised QCD - the so-called static potential 
- is simply irrelevant to the question of confinement in the real world, in which light 
quarks are ubiquitous. In fact, it is a basic feature of QCD that light-particle creation and 
annihilation effects are essentially nonperturbative and therefore it is impossible in principle 



to compute a potential between two light quarks [37.l3£ 

Drawing on a long list of sources; e.g., Refs. 39l442l|. a perspective on confinement 
was laid out in Ref. [43] ■ Confinement can be related to the analytic properties of QCD's 
Schwinger functions, which are often called Euclidean-space Green functions. For example, 



it can be read from the reconstruction theorem [4J, |45| that the only Schwinger func- 
tions which can be associated with expectation values in the Hilbert space of observables; 
namely, the set of measurable expectation values, are those that satisfy the axiom of re- 
flection positivity. This is an extremely tight constraint. It can be shown to require as a 
necessary condition that the Fourier transform of the momentum-space Schwinger function 
is a positive-definite function of its arguments. This condition suggests a practical confine- 
ment test, which can be used with numerical solutions of the DSEs (see, e.g.. Sec. III. C 



of Ref. 46|] and Sec. IV of Ref. [43]). The implications and use of reflection positivity are 
discussed and illustrated in Sec. 2 of Ref. [&\. 

It is noteworthy that any 2-point Schwinger function with an inflexion point at p^ > 
must breach the axiom of reflection positivity, so that a violation of positivity can be de- 
termined by inspection of the pointwise behaviour of the Schwinger function in momentum 
space (Sec. IV.B of Ref. fr?]). 

Consider then A{k'^), which is the single scalar function that describes the dressing 
of a Landau-gauge gluon propagator. A large body of work has focused on exposing the 
behaviour of A(A;^) in the pure Yang-Mills sector of QCD. These studies are reviewed in 



Ref. [48(. A connection with the expression and nature of conflnement in the Yang-Mills 
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IR-massive but UV-ma^ess, confined gluon 
Qgrturbative, mmslg^gluon 




massive, unconfined gluon 



k^ (GeV^) 

FIG. 2: A(fc^), the function that describes dressing of a Landau-gauge gluon propagator, plotted for 
three distinct cases. A bare gluon is described by A(fc^) = l//c^ (the dashed line), which is plainly 
convex on k^ G (0, oo). Such a propagator has a representation in terms of a non-negative spectral 
density. In some theories, interactions generate a mass in the transverse part of the gauge-boson 
propagator, so that A(fc^) — l/{k'^ + rrii), which can also be represented in terms of a non-negative 
spectral density. In QCD, however, self-interactions generate a momentum-dependent mass for the 
gluon, which is large at infrared momenta but vanishes in the ultraviolet [48| . This is illustrated by 
the curve labelled "IR-massive but UV-massless." With the generation of a mass- function, A(fc^) 
exhibits an inflexion point and hence cannot be expressed in terms of a non- negative spectral density. 



sector is indicated in Fig.[2l The appearance of an inflexion point in the two-point function 
generated by the gluon's momentum-dependent mass-function is impossible to overlook. 
Hence this gluon cannot appear in the Hilbert space of observable states. The inflexion 
point possessed by Af(p^), visible in Fig.[Tl conveys the same properties on the dressed- 
quark propagator. 

Numerical simulations of lattice-QCD confirm the appearance of an inflexion point in 
both the dressed-gluon and -quark propagators; e.g., see Fig. [1] and Ref. J48|. The signal is 
clearest for the gluon owing to the greater simplicity of simulations in the pure Yang- Mills 
We emphasise that this sense of confinement is essentially quantum field 



49|-l5l[ 



sector 

theoretical in nature. Amongst its many consequences is that confinement in QCD cannot 

veraciously be represented in potential models. 

From the perspective that confinement can be related to the analytic properties of 
QCD's Schwinger functions, the question of light-quark confinement can be translated into 
the challenge of charting the infrared behavior of QCD's universal /3-function. (Although 
this function may depend on the scheme chosen to renormalise the theory, it is unique 
within a given scheme [5^. Of course, the behaviour of the /3-function on the perturbative 
domain is well known.) This is a well-posed problem whose solution is an elemental goal of 
modern hadron physics and which can be addressed in any framework enabling the nonper- 
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turbative evaluation of renormalisation constants. It is the /3-function that is responsible 
for the behaviour evident in Figs.[T]and[2l and one of the more interesting of contemporary 
questions is whether it is possible to reconstruct the /3-function, or at least constrain it 
tightly, given empirical information on the gluon and quark mass functions. 



IV. GAP AND BETHE-SALPETER EQUATIONS 

In order to proceed it is necessary to describe explicitly the best known and sim- 
plest DSE. The Dyson or gap equation describes how quark propagation is influenced by 
interactions; viz., for a quark of flavour /, 

Sf{p)-' = Z2{h-p + m)n + Z^j g^D^^ip - q)^j^Sf{q)^Tl{q,p), (1) 

where: D^^ is the gluon propagator; Tl, the quark-gluon vertex; f , a symbol that repre- 
sents a Poincare invariant regularization of the four-dimensional Euclidean integral, with 
A the regularization mass-scale; rrif^A), the current-quark bare mass; and Zi_2(C^)^'^)) 
respectively, the vertex and quark wave function renormalisation constants, with (" the 
renormalisation point - dependence upon which is not usually made explicit. 
The gap equation's solution is the dressed-quark propagator, 

SM = 1 = ^(^^C^) ^2) 

^^' ij-pA{p-^,e) + Bip'^,e) h-p + Mip"^)' ^^ 

which is obtained from Eq. ([1]) augmented by a renormalisation condition. A mass- 
independent scheme is a useful choice and can be implemented by fixing all renormalisation 



constants in the chiral limit. (See, e.g., Ref. 53|] and references therein; or Ref. 5j] for a 
detailed discussion of renormalisation.) 

The mass function, M{p'^) = B{p'^,C,'^)/A{p'^,Q'^), is independent of the renormalisa- 
tion point, Q] and the renormalised current-quark mass. 



m 



Z^(C,A)mb-(A) = Z^^Z^mY. (3) 



wherein Z^ is the renormalisation constant associated with the Lagrangian's mass-term. 
Like the running coupling constant, this "running mass" is familiar from textbooks. How- 
ever, it is not commonly appreciated that mf' is simply the dressed-quark mass function 
evaluated at one particular deep spacelike point; viz, 

m^ = Mf{e) . (4) 

The renormalisation-group invariant current-quark mass may be inferred via 



rhf = lim 



>oo 



'Xn"" 



2 A2 



7m 

2 



Mf{p'), (5) 
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where 7^ = 12/(33 — '^Nf): Nf is the number of quark flavours employed in computing 
the running coupling; and Aqcd is QCD's dynamically-generated renormalisation-group- 
invariant mass-scale. The chiral limit is expressed by 

m/ = . (6) 



Moreover, 






VC > Aqcd, 


rrif 

h 

rrif 


_ rhf, 
rhf^ 



(7) 

However, we would like to emphasise that in the presence of DCSB the ratio mr^ jrfvr^ = 
Mf-^{p'^)/Mf^{p'^) is not independent of p^: in the infrared; i.e., Vp^ < ^qcD' i* then 
expresses a ratio of constituent-like quark masses, which, for light quarks, are two order s- 



of-magnitude larger than their current-masses and nonlinearly related to them |l5l . |l| 

The gap equation illustrates the features and flaws of each DSE. It is a nonlinear 
integral equation for the dressed-quark propagator and hence can yield much-needed non- 
perturbative information. However, the kernel involves the two-point function Df^i, and 
the three-point function Fl. The gap equation is therefore coupled to the DSEs satis- 
fied by these functions, which in turn involve higher n-point functions. Hence the DSEs 
are a tower of coupled integral equations, with a tractable problem obtained only once a 
truncation scheme is specified. It is unsurprising that the best known truncation scheme 
is the weak coupling expansion, which reproduces every diagram in perturbation theory. 
This scheme is systematic and valuable in the analysis of large momentum transfer phe- 
nomena because QCD is asymptotically free but it precludes any possibility of obtaining 
nonperturbative information. 

Given the importance of DCSB in QCD, it is significant that the dressed-quark propa- 
gator features in the axial- vector Ward-Takahashi identity, which expresses chiral symmetry 
and its breaking pattern: 

P^rg(A;; P) + I K(C) + m,(C)] vi^{k- P) = Sj\k+)ij, + ^75^;^^:-) , (8) 

where P = pi + p2 is the total-momentum entering the vertex and k is the relative- 
momentum between the amputated quark legs. To be explicit, k = {1 — rj)pi + rjp2-, with 
ry G [0, 1], and hence fe+ = pi = A; + ryP, k^ = p2 = k — [1 — r])P. In a Poincare covariant 
approach, such as presented by a proper use of DSEs, no observable can depend on r/; i.e., 
the definition of the relative momentum. NB. Sec . IV V^ discusses the important differences 
encountered in treating fiavourless pseudoscalar mesons. 

In Eq. ([8]), Vl^f^ and Fg are, respectively, the amputated axial- vector and pseudoscalar 
vertices. They are both obtained from an inhomogeneous Bethe-Salpeter equation (BSE), 
which is exemplified here using a textbook expression [55(: 

[r5;.(A:; P)]tu = Z2h5%.]tu + [ [Siq+)r,^iq; P)5(g_)],,A'[„^(<?, k; P), (9) 

Ja 
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in which K is the fully-amputated quark-antiquark scattering kernel, and the colour-, Dirac- 
and flavour-matrix structure of the elements in the equation is denoted by the indices 
r, s, t, u. NB. By definition, K does not contain quark-antiquark to single gauge-boson 
annihilation diagrams, nor diagrams that become disconnected by cutting one quark and 
one antiquark line. 

The Ward-Takahashi identity, Eq. ^, entails that an intimate relation exists between 
the kernel in the gap equation and that in the BSE. (This is another example of the coupling 
between DSEs.) Therefore an understanding of chiral symmetry and its dynamical breaking 
can only be obtained with a truncation scheme that preserves this relation, and hence 
guarantees Eq. ([8]) without a fine-tuning of model-dependent parameters. 



V. NONPERTURBATIVE TRUNCATION 

Through the gap and Bethe-Salpeter equations the pointwise behaviour of the /3- 
function determines the pattern of chiral symmetry breaking; e.g., the behaviour in Fig.[TJ 
Moreover, the fact that these and other DSEs connect the /3- function to experimental 
observables entails that comparison between computations and observations of the hadron 
mass spectrum, and hadron elastic and transition form factors, can be used to constrain 
the /3-function's long-range behaviour. 

In order to realise this goal, a nonperturbative symmetry-preserving DSE truncation 
is necessary. Steady quantitative progress can be made with a scheme that is systematically 



improvable |56l . l57l |. In fact, the mere existence of such a scheme has enabled the proof of 
exact nonperturbative results in QCD. 

Before describing a number of these in some detail, it is worth explicating the range 



of applications. For example, there are: veracious statements about the pion a-terni 15l |: 
radially-excited and hybrid pseudoscalar mesons 5a, |59(] ; heavy-light [60] and heavy-heavy 
mesons [30y ; novel results for the pion susceptibility obtained via analysis of the isovector- 



pseudoscalar vacuum polarisation 6l[ , which bear upon the essential content of the so-called 
"Mexican hat" potential that is used in building models for QCD; and a derivation [621] of 
the Weinberg sum rule [6;^. 

V-1. Pseudoscalar meson mass formula 

Turning now to a fuller illustration, the first of the results was introduced in Ref. 64i |: 
namely, a mass formula that is exact for flavour non-diagonal pseudoscalar mesons: 

Ih^- mjf^_ = (mj^ + mj^)pii^_ > (10) 

where: m} are the current-masses of the quarks constituting the meson; and 

Ji 

fHg_P^^ = {0\qf2-f5-/f,qf^\Ho~^) = ZatrcD/ i-f5ltMSf^{q+)TH^_{q; P)Sf2{q^) , (11) 

Jq 

ipHg- = -{0\qf2il5qh\Ho-) = ZitrcD 75%(g+)r//^,_ (g; P)5'/2(g_) , (12) 

Ja 
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where T w _ is the pseudoscalar meson's bound-state Bethe-Salpeter ampHtude: 

Th,. {k; P) = 75 [iEH,^ {k; P) + j ■ PFh^^_ {k; P) 

+7 • k Gh^^ (k; P) - a^.k^P.Hn^^^ {k; P)] , (13) 

which is determined from the associated homogeneous BSE. 

It is worth emphasising that the quark wavefunction and Lagrangian mass renormal- 
isation constants, Z2,4(C,A), respectively, depend on the gauge parameter in precisely the 
manner needed to ensure that the right-hand sides of Eqs. (jlip . (J12p are gauge- invariant. 
Moreover, Z2(C, A) ensures that the right-hand side of Eq. ([TT]) is independent of both C 
and A, so that fn _ is truly an observable; and Z4(C, A) ensures that p^ is independent 

of A and evolves with ^ in just the way necessary to guarantee that the product m^p^^ is 
renormalisation-point-independent. In addition, it should be noted that Eq. (|10|) is valid for 
every pseudoscalar meson and for any value of the current-quark masses; viz., rhf^ £ [0, oo), 
i = 1,2. This includes arbitrarily large values and also the chiral limit, in whose neigh- 



bourhood Eq. (jlOp can be shown 6J] to reproduce the familiar Gell-Mann-Oakes-Renner 
relation. 

The axial- vector Ward-Takahashi identity, Eq. (jH]), is a crucial bridge to Eqs. pO|) - 
(I12p : and on the way one can also prove the following Goldberger-Treiman-like relations 

FR{k;0) + 2fH^_FH^_{k;0) = A''{k^), (15) 

Gnik; 0) + 2ffj^_ Gh^_ (k; 0) = ^^°(fe') , (16) 

HR{k;0)+2fH^^_HH^_{k;0) = 0, (17) 

wherein the superscript indicates that the associated quantity is evaluated in the chiral 
limit, and Fn, Gr, Hr are analogues in the inhomogeneous axial- vector vertex of the scalar 
functions in the Hq— meson's Bethe-Salpeter amplitude. 

These identities are of critical importance in QCD. The first, Eq. (|14p . can be used to 
prove that a massless pseudoscalar meson appears in the chiral-limit spectrum if, and only 
if, chiral symmetry is dynamically broken. Moreover, it exposes the fascinating consequence 
that the solution of the two-body pseudoscalar bound-state problem is almost completely 
known once the one-body problem is solved for the dressed-quark propagator, with the 
relative momentum within the bound-state identified unambiguously with the momentum 
of the dressed-quark. This latter emphasises that Goldstone's theorem has a pointwise 
expression in QCD. The remaining three identities are also important because they show 
that a pseudoscalar meson must contain components of pseudovector origin. This result 
overturned a misapprehension of twenty-years standing; namely, that only Eh _ (A;; 0) is 
nonzero |65J]. These pseudovector components materially influence the observable properties 
of pseudoscalar mesons GGJ, I68l47ll] , as do their analogues in other mesons 72h74|| . 
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FIG. 3: Lattice-QCD results for the ratio of the decay constants for the first-excited- and ground- 
state pseudoscalar mesons as a function of the pion mass squared. (Lattice parameters: volume = 
16"^ X 32; /3 = 5.2, spacing a ~ 0.1 fm, two flavours of degenerate sea quarks; Wilson gauge action and 
clover fermions.) The "not improved" results were obtained from a fermion action with poor chiral 
symmetry properties. In this case I/tti/ZttoI ~ 0-4, consistent with expectations based on quantum 
mechanics. The "improved" results were obtained through implementation of the full ALPHA 
method for the nonperturbative improvement of the fermion action, which greatly improves the 
simulation's chiral symmetry properties. In this case, IfTTi/f-Kol ~ 0.01. (NB. The sign of the ratio 
was not determined in the simulation but is discussed in Ref. i77j. Figure adapted from Ref. [78].) 



It is natural to reiterate here a prediction for the properties of non-ground-state 
pseudoscalar mesons, which fohows from the exact results described above; namely, in the 
chiral limit 5a, l59[] 



A„=0,Vn>l. 



(18) 



This is the statement that Goldstone modes are the only pseudoscalar mesons to possess 
a nonzero leptonic decay constant in the chiral limit when chiral symmetry is dynamically 
broken. The decay constants of all other pseudoscalar mesons on this trajectory, e.g., radial 
excitations, vanish. On the other hand, in the absence of DCSB the leptonic decay constant 
of each such pseudoscalar meson vanishes in the chiral limit; i.e, Eq. (llSp is true Vn > 0. 

From the perspective of quantum mechanics, Eq. (jlSp is a surprising fact. The lep- 
tonic decay constant for 5-wave states is typically proportional to the wave function at the 
origin. Compared with the ground state, this is smaller for an excited state because the wave 
function is broader in configuration space and wave functions are normalised. However, it 
is a modest effect; e.g., consider the e~^e~ decay of vector mesons, for which a calculation 
in relativistic quantum mechanics based on light-front dynamics [75( yields \fpi/fpo\ = 0.5, 



consistent with the value inferred from experiment. Thus, it is not uncommon for Eq. ([18 
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to be perceived as "remarkable" or "unbelievable." Notwithstanding this, in connection 



with the pion's first radial excitation, the value of /t^ = — 2MeV predicted in Ref. 58|] is 
consistent with experiment \76\ and simulations of lattice-QCD [78i], as illustrated in Fig. [31 
It is now recognised that the suppression of /jr^ is a useful benchmark, which can be used 
to tune and validate lattice QCD techniques that try to determine the properties of excited 
states mesons. 

V-2. In-hadron condensates 

For the last thirty years, condensates; i.e., nonzero vacuum expectation values of local 
operators, have been used as parameters in order to correlate and estimate essentially non- 
perturbative strong-interaction matrix elements [79|]. They are also basic to current algebra 
analyses. It is conventionally held that such quark and gluon condensates have a physical 
existence, which is independent of the hadrons that express QCD's asymptotically realis- 
able degrees-of-freedom; namely, that these condensates are not merely mass-dimensioned 
parameters in a theoretical truncation scheme, but in fact describe measurable spacetime- 
independent configurations of QCD's elementary degrees-of-freedom in a hadronless ground 



state. Owing to confinement, however, this view is erroneous [34i |. 

One may readily recapitulate upon the argument. To begin, note that Eq. (jllh is the 
exact expression in QCD for the leptonic decay constant of a pseudoscalar meson. It is 
a property of the pion and, as consideration of the integral expression reveals, it can be 
described as the pseudovector projection of the pion's Bethe-Salpeter wavefunction onto 
the origin in configuration space. Note that the product ip = STS is called the Bethe- 
Salpeter wavefunction because, when a nonrelativistic limit can validly be performed, the 
quantity V' at fixed time becomes the quantum mechanical wavefunction for the system 
under consideration. (NB. In the neighborhood of the chiral limit, a value for fn _ can be 



estimated via either of two approximation formulae [SOHSj]. These formulae both illustrate 
and emphasize the role of fn _ as an order parameter for DCSB.) 

If chiral symmetry were not dynamically broken, then in the neighborhood of the 
chiral limit /// _ oc m |58t]. Of course, chiral symmetry is dynamically broken in QCD 



12h14|] and for the ground-state pseudoscalar 



lim/^^.(m) = /0 ^0. (19) 

Taken together, these last two observations express the fact that fn _ , which is an intrinsic 
property of the pseudoscalar meson, is a bona fide order parameter for DCSB. An analysis 
within chiral perturbation theory [83| suggests that the chiral limit value, Z]^ _ , is ~ 5% 
below the measured value of 92.4 MeV; and efficacious DSE studies give a 3% chiral-limit 



reduction [Sjj. 



Now, Eq. ()12p is kindred to Eq. (jlip : it is the expression in quantum field theory 
which describes the pseudoscalar projection of the pseudoscalar meson's Bethe-Salpeter 
wavefunction onto the origin in configuration space. It is thus truly just another type of 
pseudoscalar meson decay constant. 

In this connection it is therefore notable that one may rigorously define an "in-meson" 
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condensate; viz. J64l. |84| : 



-.J 



k)j (m) . 



Now, using Eq. ^^, one finds [6j] 



lim 4 (m) = Z4trcD / 44-,S\q;0 = -{qq)l 



(20) 



(21) 



Hence the so-called vacuum quark condensate is, in fact, the chiral-hmit value of the in- 
meson condensate; i.e., it describes a property of the chiral-limit pseudoscalar meson. One 
can therefore argue that this condensate is no more a property of the "vacuum" than the 
pseudoscalar meson's chiral-limit leptonic decay constant. Moreover, Ref. |;85] establishes 
the equivalence of all three definitions of the so-called vacuum quark condensate: a constant 
in the operator product expansion [8£ 



871 ]: via the Banks-Casher formula 88|; and the 



trace of the chiral-limit dressed-quark propagator. (NB. In the presence of confinement it is 
impossible to write a valid nonperturbative definition of a single quark or gluon annihilation 
operator and therefore impossible to rigorously define a second quantised vacuum for QCD. 
To do so would be to answer the question: What is the state that is annihilated by an 
operator which is unknowable?) 

The chiral-limit vacuum quark condensate is therefore qualitatively equivalent to the 
pseudoscalar-meson leptonic decay constant in the sense that both are obtained as the 
chiral-limit value of well-defined gauge-invariant hadron-to-vacuum transition amplitudes 
that possess a spectral representation in terms of the current-quark- mass. Thus, whereas it 
might sometimes be convenient to imagine otherwise, neither is essentially a constant mass- 
scale that fills all spacetime. Hence, in particular, the quark condensate can be understood 
as a property of hadrons themselves - a perspective also advocated in Ref. [89[ and now 
established for all hadrons [36[] - which is expressed, for example, in their Bethe-Salpeter 
or light-front wave functions. In the latter instance, the light-front-instantaneous quark 
propagator appears to play a crucial role [34, |90(] . 

This has enormous implications for the cosmological constant. The universe is ex- 
panding at an ever-increasing rate and theoretical physics has tried to explain this in terms 
of the energy of quantum processes in vacuum; viz., condensates carry energy and so, if they 
exist, must contribute to the universe's dark energy, which is expressed in the computed 
value of the cosmological constant. The problem is that hitherto all potential sources of 
this so-called vacuum energy give magnitudes that far exceed the value of the cosmological 
constant that is empirically determined. This has been described as "the biggest embar- 



rassment in theoretical physics" 9l[. However, given that, in the presence of confinement. 



condensates do not leak from within hadrons, then there are no space-time-independent 
condensates permeating the universe 3J, |3g] . This nullifies completely their contribution 
to the cosmological constant and reduces the mismatch between theory and observation by 
a factor of 10^^ [35| , and possibly by far more, if technicolour-like theories are the correct 
paradigm for extending the Standard Model. 
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FIG. 4: This simple flavour-mixing diagram, is immaterial to the resolution of the 77-77' conundrum, 
as is any collection of processes for which the figure may serve as a skeleton diagram. (Straight lines 
denote quarks and springs denote gluons.) 



V-3. Flavourless pseudoscalar mesons 

In connection with electric-charge-neutral pseudoscalar mesons, Eq. (|10p is strongly 
modified owing to the non-Abelian anomaly. This entails that whilst the classical action as- 
sociated with QCD is invariant under UA{Nf) (non-Abelian axial transformations generated 
by Ao75, where Aq oc diag[l, . . . , Iat^]), the quantum field theory is not. The modification 
is particularly important to properties of i] and r]' mesons. The latter is obviously a pecu- 
liar pseudoscalar meson because its mass is far greater than that of any other light-quark 
pseudoscalar meson; e.g., m^/ = 1.75 m^. We note that the diagram depicted in Fig. H] is 
often cited as central to a solution of the 77-77' puzzle. However, as will become clear below, 
whilst it does contribute to flavour-mixing, the process is immaterial in resolving the 77-77' 
conundrum. 

The correct mass formula for flavourless pseudoscalars follows from consideration of 
the complete C/a(-^/) Ward-Takahashi identity: 

P^j:l^{k- P) = S'Hk+)z^,T'' + i-f,T-S-Hk^) - 2iM'''rl{k; P) - A^{ 



.P) 



(22) 



which generalises Eq. ([8]). In Eq. (j22p . {J-'^\ a = 0, . . . , N^ — 1} are the generators of U{Nf) in 
the fundamental representation, orthonormalised according to trJ^°J^'' = 2"'i the dressed- 
quark propagator S = diag[S'u, Sd, Ss,Sc, S^,, ■ ■ ■] is matrix- valued; and 



M 



ab 



tVp 



{-F^A^^}^* 



(23) 



where Ai'^ is a matrix of renormalised current-quark masses and the trace is over flavour 
indices. 

The final term in the last line of Eq. (|22p expresses the non-Abelian axial anomaly. 
It can be written 



A''ik;P) 
Au{k;P) 



S-\k+)6''°Au{k;P)S-\k^), 

fd^xd^y e'^^+'^'-^-'y^Nf {J^ q{x) Q(0) q{y)) , 

and since A'^^^{k; P) is a pseudoscalar, it has the general form 

A\k-P) = ^^5[i£A{k;P) 

+7 • PJ'Ai.k; P)+j-kk- PGAik; P) + a^,k^P,nA{k; P)] 



(24) 
(25) 



(26) 
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The matrix element in Eq. (j25p represents an operator expectation value in full QCD; the 
operation in Eq. (|24p amputates the external quark lines; and 



ry 

Q{x) = i-^tic [e^upaF^uFpaix)] = d^K^ix) (27) 

47r 

is the topological charge density operator, where the trace is over colour indices and Ffj_^ = 
2X"'F!^^ is the matrix-valued gluon field strength tensor. It is plain and important that 
only A^^^ is nonzero. NB. While Q{x) is gauge invariant, the associated Chern-Simons 
current, K^, is not. Thus in QCD no physical state can couple to K^ and hence no state 
which appears in the observable spectrum can contribute to a resolution of the so-called 
C/A(l)-pi'oblem; i.e., physical states cannot play any role in ensuring that the r]' is not a 
Goldstone mode. 

As described in Sec lVVH^ if one imagines there are Nf massless quarks, then DCSB 
is a necessary and sufficient condition for the a 7^ components of Eq. (I22p to guarantee 
the existence of N'i — 1 massless bound-states of a dressed-quark and -antiquark. However, 
owing to Eq. (I24p . a = in Eq. (j22p requires special consideration. One case is easily 
covered; viz., it is clear that if A^ = 0, then the a = component of Eq. (j22p is no different 
to the others and there is an additional massless bound-state in the chiral limit. 

On the other hand, the large disparity between the mass of the r/'-meson and the 
octet pseudoscalars suggests that ^'^ 7^ in real- world QCD. If one carefully considers that 



possibility, then the Goldberger-Treiman relations in Eqs. (|14p - (|17p become 93] 

'ifH,.EBs{k;0) = 2B\k^)-£,^{k;0), (28) 

F^ik;0) + 2fH,.FBsik;0) = A^ik^) - :F_^{k;0), (29) 

G%ik;0) + 2ffj^^_GBsik;0) = 2 ^^"(A:^) - g^(fc;0), (30) 

H%{k;0) + 2fH^^_HBs{k;0) = -?^^(A;;0), (31) 

It follows that the relationship 

SAk;0)=2B'ik^), (32) 

is necessary and sufficient to guarantee that Eg (A;;P), the flavourless pseudoscalar vertex, 
does not possess a massless pole in the chiral limit; i.e., that there are only N'j — 1 massless 
Goldstone bosons. Now, in the chiral limit, B^(k'^) 7^ if, and only if, chiral symmetry 
is dynamically broken. Hence, the absence of an additional massless bound-state is only 
assured through the existence of an intimate connection between DCSB and an expectation 
value involving the topological charge density. 

This critical connection is further highlighted by the following result, obtained 
through a few straightforward manipulations of Eqs. ([2^ . ([M|) and ([2^ : 

{qq)l = - hm 4 (m) = - lim Z4(C', A^) trcD / S^{q,C) (33) 

= ^ I d^x{q{x)i^,q{x)Qm^. (34) 
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The absence of a Goldstone boson in the a = channel is only guaranteed if this explicit 
identity between the chiral-limit in-nieson condensate and a mixed vacuum polarisation 
involving the topological charge density is satisfied. 



Mass formulae valid for all pseudoscalar mesons have also been obtained [94] 
Ih,. rn\_ = 2 M^'p'h^^ + 6^' nH^_ , (35) 

where 

/^^_ P^ = Z2tr f ^'^tsTm Xh„_ ((?; P) , (36) 

Jq 

ip%^^_{0 = Z^tr f J^^J^Xh.Aq^P), (37) 

Jq 

UH,, = J^^H,^ , m,. = {m\Ho-) ■ (38) 



For charged pseudoscalar mesons, Eq. (p5|) is equivalent to Eq. pOj) . but the novelty of 
Eq. ()35p is what it expresses for flavourless pseudoscalar s. To illustrate, consider the case 
of a U{Nf = 3)-symmetric mass matrix, in which all Nf = 3 current-quark masses assume 
the single value m^, then this formula yields: 

Plainly, the t]' is split from the Goldstone modes so long as n^i 7^ 0. Numerical simulations 
of lattice-QCD have confirmed this identity 93|, l94| - 



It is important to elucidate the physical content of n,,'. Returning to the definition: 

liolQW), (40) 

it is readily seen to be another type of in-meson condensate. It is analogous to those 
discussed in Sec. I VV^ but in this case the hadron-to- vacuum transition amplitude measures 
the topological content of the r]' . One may therefore state that the r/' is split from the 
Goldstone modes so long as its wavefunction possesses nonzero topological content. This is 
plainly very different to requiring that the QCD vacuum is topologically nontrivial. 

Within QCD the properties of the r]' can be computed via the BSE, just like other 
mesons. A nonzero value of rirj' can be achieved with a Bethe-Salpeter kernel that contains 
the contribution depicted in Fig.O because one may argue from Eqs. ([2S]) and ([27]) that an 
anomaly-related contribution to a meson's Bethe-Salpeter kernel cannot contain external 
quark or antiquark lines that are connected to the incoming lines: purely gluonic configu- 
rations must mediate, as illustrated in Fig.[5j Furthermore, it is straightforward to see that 
no finite sum of gluon exchanges can serve this purpose. Indeed, consider any one such 
single contribution in the chiral limit. It will be proportional to the total momentum and 
hence vanish for P = 0, in conflict with Eq. (I39p . This lies behind the need for something 
like the Kogut-Susskind ghost; i.e., the coupling of a massless axial- vector gauge-like field to 
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FIG. 5: An illustration of the nature of the contribution to the Bethe-Salpeter kernel associated 
with the non-Abelian anomaly. All terms have the "hairpin" structure illustrated in the lower panel. 
No finite sum of such intermediate states is sufhcient. (Straight lines denote quarks, with /i and ji 
independent, and springs denote gluons.) (Figure adapted from Ref. 92].) 



the Chern-Simons current, which does not ap pea r in the particle spectrum of QCD because 
the current is not gauge invariant. (See Ref. 95|] and Sec. 5.1 of Ref. J96i].) 
It is argued [97,, >M] that in QCD with N^ colours, 



n^, ~ 



^Wr 



(41) 



and it can be seen to follow from the gap equation, the homogeneous BSE and Eqs. (j36p . 
(1371) that 



Jr]' 



p'AC). 



One thus obtains 



n„i 



m^ 



In' 



fiAC) 






(42) 



(43) 



The first term vanishes in the limit Nc ^ oo while the second remains finite. Subsequently 
taking the chiral limit, the r]' mass approaches zero in the manner characteristic of all 
Goldstone modes. (NB. One must take the limit Nc ^ oo before the chiral limit because 
the pro cedu res do not commute {99f|.) These results are realised in the effective Lagrangian 
of Ref. [100] in a fashion that is consistent with all the constraints of t he anoma lous Ward 
identity. NB. This is not true of the so-called 'tHooft determinant 96l. llOL 102l |. 

Implications of the mass formula in Eq. (j35p were exemplified in Ref. [9j] using an ele- 
mentary dynamical model that includes a one-parameter Ansatz for that part of the Bethe- 
Salpeter kernel related to the non-Abelian anomaly, an illustration of which is provided in 
Fig.m The study compares ground-state pseudoscalar- and vector-mesons constituted from 
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all known quarks, excluding the t-quark. Amongst the notable results is a prediction for 
the mixing angles between neutral mesons; e.g., 

On = -15.4° , On' = -15.7° . (44) 

NB. There are necessarily two mixing angles, with each determined at the appropriate pole 
position in the inhomogeneous vertex. It is interesting that the angles are approxim ately 



equal and compare well with the value inferred from a single mixing angle analysis 103l | 
9 = -13.3° ±1.0°. 

It is worth explicating the nature of the flavour-induced difference between the vr^ 
and vr masses. If one ignores mixing with mesons containing other than u, d-quarks; viz., 
works solely within SU{Nf = 2), then m^o — m^+ = — 0.04MeV. On the other hand, the 
full calculation yields m^o — m^+ = — 0.4MeV, a factor of ten greater, and one obtains a 
TT^-T] mixing angle, whose value at the neutral pion mass shell is 



'Tvr] 



[ml.) = 1-2°. (45) 



For comparison, Ref. [104i | infers a mixing angle of 0.6° it 0.3° from a ii'-matrix analysis 
of the process pd — )• ^Hevr''. Plainly, mixing with the r/- meson is the dominant non- 
electromagnetic contribution to the vr -vr'' mass splitting. The analogous angle at the rj 
mass-shell is 

0^^(m^) = 1.3°. (46) 

The angles in Eq. (|44p correspond to 

|7r°) ~ 0.72 uu- 0.69 Jd- 0.013 ss, (47) 

\r]) ~ 0.53UU + 0.57 dd-Omss, (48) 

\rj') ~ 0.44 tZu + 0.45 fi(i + 0.78 ss. (49) 

Evidently, in the presence of a sensible amount of isospin breaking, the vr'^ is still predom- 
inantly characterised by J-"'^ but there is a small admixture of ss. It is found in Ref. [92] 
that mixing with the vr^ has a similarly modest impact on the flavour content of the 7] and 
T]' . It's effect on their masses is far less. 

VI. EXPRESSING DCSB IN BOUND-STATES 

Despite the successes achieved with the systematic scheme, its practical application 
has numerous shortcomings. For example, the leading-order truncation is accurate for 
electrically-charged ground-state pseudoscalar- and vector-mesons because corrections in 
these channels largely cancel, owing to parameter-free pr eserv ation of the Ward-Takahashi 



identities. However, they do not cancel in other channels [105l4l08] . Hence studies based on 
the rainbow-ladder truncation, or low-order improv ements th ereof, have usually provided 
poor results for scalar- and axial- vector- mesons 321. ll09Mll3l |. produced masses for exotic 
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states that are too low in comparison with other estimates 32, 1771 . 1 103 . Ill3 |. and e xhibi t 
gross sensitivity to model parameters for tensor-mesons 1141 ] and excited states 58l. l77l . lll5l |. 
In these circumstances one must conclude that physics important to these states is omitted. 
One anticipates therefore that significant qualitative advances in understanding the essence 
of QCD could be made with symmetry-preserving kernel Ansdtze that express important 
additional nonperturbative effects, which are impossible to capture in any finite sum of 
contributions. Such an approach has recently become available 
summarising herein. 



82( 1 and is well worth 



VI-1. Building the Bethe-Salpeter kernel 

Consider, e.g., flavoured pseudoscalar and axial- vector mesons, which appear as poles 
in the inhomogeneous BSE for the axial-vector vertex, F^, where f,g are flavour labels. 
An exact form of that equation is {k, q are relative momenta, P is the total momentum 
flowing into the vertex, and q± = q ± P/2, etc.) 

^f9^u.m - ^— f'^9'D^p{k-q)^^^Sf{q+)rlUq-,P)S,{q^)^rUq^ 



Vmk-P) = Z2757M 



_,A:_ 



+ 



g^D, 



A'^ 



A" 



{k-q)-^^Sf{q+)-4Uk,q;P), 



(50) 



where A^^o is a 4-point Schwinger function. (The pseudoscalar verte x satisfies an analogue 



of Eq. (I50p .) This form of the BSE was first written in Ref. |l07l | and is illustrated in 
the lower-panel of Fig. ([6]) . The diagrammatic content of the right-hand-side is completely 
equivalent to that of Eq. ([9]), which is depicted in the upper-panel of the figure. However, 
in striking qualitative opposition to that textbook equation, Eq. ([50|) partly embeds the 
solution vertex in the four-point function, A, whilst simultaneously explicating a part of 
the effect of the dressed-quark-gluon vertex. This has the invaluable consequence of enabling 
the derivation of both an integral equation for the n ew B ethe-Salpeter kernel. A, in which 
the driving term is the dressed-quark-gluon vertex 1071 ]. and a Ward-Takahashi identity 
relating A to that vertex [82]. No similar equations have yet been found for K and hence 
the textbook form of the BSE, whilst tidy, is very limited in its capacity to expose the 
effects of DCSB in hadron physics. 

As emphasised above, no study of light-quark hadrons is dependable if it fails to 
comply with the axial- vector Ward-Takahashi identity, Eq. ([8]). The condition 

P^Ag^(A:,Q;P) + i[mj(C) + m,(C)]Ag(A:,(z;P)=F^((?+,A:+)z75 + i75r^(g-,fc-), (51) 

where A^? is the analogue of A^^o in the pseudoscalar equation, is necessary and sufficient 
to ensure the Ward-Takahashi identity is satisfied by the solution of Eqs. ([1]) and (j50|) [821]. 
Consider Eq. (I5ip . Rainbow-ladder is the leading-order term in the systematic DSE 



truncation scheme of Refs. [56|, |57l]. It corresponds to Ti = ^y, in which case Eq. ([5T]) is 
solved by A^^o = = Ag?. This is the solution that indeed provides the rainbow-ladder 
forms of Eq. (j50p . Such consistency will be apparent in any valid systematic term-by-term 
improvement of the rainbow-ladder truncation. 
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FIG. 6: Omitting the inhomogeneity, the upper panel illustrates the textbook form of the Bethe- 
Salpeter equation, Eq. ([9]), whereas the lower panel depicts the form expressed in Eq. (l50t . The 
reversal of the total-momentum's flow is immaterial here. NB. In any symmetry-preserving trunca- 
tion, beyond the leading-order identified in Ref. [57|; i.e., r ainbo w-ladder, the Bethe-Salpeter ke rnel is 
nonplanar even if the vertex in the gap equation is planar [lOTJ. (Figure adapted from Refs. |ll. ll07l |.) 



However, Eq. (j5ip is far more than just a device for checking a truncation's consis- 
tency. For, just as the vector Ward-Tak ahashi id entity has long been used to build Ansdtze 
for the dressed-quark-photon vertex [l|, Ilia . Ill7l | , Eq. (i5T]l provides a tool for constructing 
a symmetry-preserving kernel of the BSE that is matched to any reasonable Ansatz for the 
dressed-quark-gluon vertex which appears in the gap equation. With this powerful capac- 
ity, Eq. (j5ip achieves a goal that has been sought ever since the Bethe-Salpeter equation 
was introduced [55|. As will become apparent, it produces a symmetry-preserving kernel 
that promises to enable the first reliable Poincare invariant calculation of the spectrum of 
mesons with masses larger than 1 GeV. 

The utility of Eq. (j5ip was illustrated in Ref. [_82;] through an application to ground 
state pseudoscalar and scalar mesons composed of equal-mass u- and d-quarks. To this end, 
it was supposed that in Eq. ([1]) one employs an Ansatz for the quark-gluon vertex which 
satisfies 



p,irak+,k^) = B{p^ 



Sj\k, 



sjHk^] 



(52) 



with B flavour-independent. (NB. While the true quark-gluon vertex does not satisfy this 
identity, owing to the form of the Slavnov- Taylor identity which it does satisfy, it is plausible 
that a solution of Eq. ()52p can provide a reasonable pointwise approximation to the true 
vertex.) Given Eq. ([52]). then Eq. ([5T]) entails {I = q — k) 



il04l{k, q; P) = B{lf Ti\q- P) - vi^k; P) 



(53) 



with an analogous equation for P^j^lpiKi^^ g{k,q]P). This identity can be solved to obtain 



Ag(A:, q- P) := B{{k - qf) 75 A^^(A:, q- P) , 



(54) 
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with, using an obvious analogue of Eq. (J13p . 

AJ,'{k,q;P)=2£p[iAE,{q,k;P)+-fPAF,{q,k;P)]+^p^Gs{q,k;P) 
+ 2e^-f£AGMk;P) + hi3n-P]^Hs{q,k;P) + 2£^[j-£,j-P]AH,{q,'<^;P), (55) 

where £ = {q + k)/2, i:^{q,k;P) = [<^{q;P) + ^{k;P)]/2 and A^{q,k;P) = [<^{q;P)- 
$(fc;P)]/[g2-A;2]. 

Now, given any Ansatz for the quark-gluon vertex that satisfies Eq. (|52p . then the 
pseudoscalar analogue of Eq. ([50|) . and Eqs. ([1]), ([M|) . (j55]) provide a symmetry-preserving 
closed system whose solution predicts the properties of pseudoscalar mesons. The relevant 
scalar meson equations are readily derive d. (NB. T he role played by resonant contributions 
to the kernel in the scalar channel la, IllSl . Ill9l | is not being overlooked but they are 



not pertinent here. Further comments appear in Sec. lVI VI-3[ ) With these systems one 
can anticipate, elucidate and understand the influence on hadron properties of the rich 
nonperturbative structure expected of the fully-dressed quark-gluon vertex in QCD: in 
particular, that of the dynamically generated dressed-quark mass function, whose impact 



is quashed at any finite order in the truncation scheme of Ref . [57 1 . 

To proceed one need only specify the gap equation's kernel because, as noted above, 
the BSEs are completely defined therefrom. To complete the illustration [8z] a simplified 
form of the effective interaction in Ref. ^] was employed and two vertex Ansdtze were 
compared; viz., the bare vertex F^ = 7^, which defines the rainbow- ladd er truncation of 
the DSEs and omits vertex dressing; and the Ball-Chiu (BC) vertex lla ]. which nonper- 
turbatively incorporates vertex dressing associated with DCSB: 

iF^((/, k) = iJ^A^ {q\ k") 7^ + 2£^ [?7 • £ Aa. {q\ k") + Ab. {q\ A:^)] . (56) 

A particular novelty of the study is that one can calculate the current-quark-mass- 
dependence of meson masses using a symmetry-preserving DSE truncation whose diagram- 
matic content is unknown. That dependence is depicted in Fig.[7] and compared with the 
rainbow-ladder result. The m-dependence of the pseudoscalar meson's mass provides nu- 
merical confirmation of the algebraic fact that the axial-vector Ward-Takahashi identity 
is preserved by both the rainbow-ladder truncation and the BC-consistent Ansatz for the 
Bethe-Salpeter kernel. The figure also shows that the axial-vector Ward-Takahashi identity 
and DCSB conspire to shie ld the pio n's mass from material variation in response to dressing 



the quark-gluon vertex M, llOTI . llOSi ] . 



As noted in Ref. |82l | , a rainbow-ladder kernel with realistic interaction strength yields 
RL._ 2M(0) - m„ 



e 



2M(0) 



= (0.3 ±0.1), (57) 

RL 



which can be contrasted with the value obtained using the BC-consistent Bethe-Salpeter 
kernel; viz.. 



-BC 



<0.1. (58) 
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FIG. 7: Dependence of pseudoscalar (left panel) and scalar (right) meson masses on the current- 
quark mass, m. The Ball-Chiu vertex (BC) result is compared with the rainbow-ladder (RL) result. 
(Figure adapted from Ref. 8^.) 



Plainly, significant additional repulsion is present in the BC-consistent truncation of the 
scalar BSE. 

Scalar mesons are commonly identified as ^Pq states. This assignment expresses a 
constituent-quark-model perspective, from which a J^^ = 0^^ fermion-antifermion bound- 
state must have the constituents' spins aligned and one unit of constituent orbital angular 
momentum. Hence a scalar is a spin and orbital excitation of a pseudoscalar meson. Of 
course, no constituent-quark-model can be connected systematically with QCD. Neverthe- 
less, as observed in Sec. [Til the presence of orbital angular momentum in a hadron's rest 
frame is a necessary consequence of Poincare covariance and the momentum-dependent 
vector-boson-exchange character of QCD [sl. ISQ. l31|. so there is a realisation in QCD of the 
quark-model anticipation. 

Extant studies of realistic corrections to the rainbow-ladder truncation show that 



they reduce hyperfine splitting 108l |. Hence, with the comparison between Eqs. (I57p and 
(I58p one has a clear indication that in a Poincare covariant treatment the BC-consistent 
truncation magnifies spin-orbit splitting, an effect which can be attributed to the influence 
of the quark's dynamically-enhanced scalar self-energy [^ in the Bethe-Salpeter kernel. 

VI-2. Quark anomalous chromomagnetic moment 

It was conjectured in Ref. [82| that a full realisation of DCSB in the Bethe-Salpeter 
kernel will have a big impact on mesons with mass greater than 1 GeV. Moreover, that it can 
overcome a longstanding failure of theoretical hadron physics. Namely, no extant continuum 
hadron spectrum calculation is believable because all symmetry preserving studies produce 
a splitting between vector and axial-vect or niesori s that is far too small: just one-quarter 
of the experimental value (see, e.g., Refs. llCll4ll2l |). Foll owing that conjecture, there have 
been significant developments, which will now be related 1231 ]. 

In Dirac's relativistic quantum mechanics, a fermion with charge q and mass m, inter- 
acting with an electromagnetic field, has a magnetic moment ^ = q/[2m]. This prediction 
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held true for the electron un til im provements in experimental techniques enabled the dis- 
covery of a small deviation 1201], with the moment increased by a multiplicative factor: 
1.00119 lb 0.00005. This correction was explained by the first systematic computation using 
renormalised quantum electrodynamics (QED) [121i] : 



q 

2m 



1 + 



a 
2^ 



Q 
2m ' 



(59) 



where a is QED's fine structure constant. The agreement with experiment established 
quantum electrodynamics as a valid tool. The correction defines the electron's anom,alous 
magnetic mom ent, which is now known with extraordinary precision and agrees with theory 
at 0{a^) [l2^ . 



The fermion-photon coupling in QED is described by: 



/ 



d^ iq 'iIj{x)j^^{x) A^ (x) 



(60) 



where ip{x), ^{x) describe the fermion field and Af^{x) describes the photon. This in- 
teraction generates the following electromagnetic current for an on-shell Dirac fermion 

{k = pf -Pi), 



iqu{pf) 



7m^i( 



^2^"''' 



KF2{k^ 



u{pi) , 



(61) 



where: -Fi(fe^), F2{k'^) are form factors; and u{p), u{p) are electron spinors. Using their 
Euclidean space definition, one can derive a Gordon-identity; viz., with 2£ = pf + pi, 



2mu{pf)ij^u{pi) = u{pf) [2i^ 



~r l<y ^yky 



u{pi 



(62) 



With this rearrangement one sees that for massive fermions the interaction can be decom- 
posed into two terms: the first describes the spin-independent part of the fermion-photon 
interaction, and is common to spin-zero and spin-half particles, whilst the second expresses 
the spin-dependent, helicity hipping part. Moreover, one reads from Eqs. (j6ip and (j62p 
that a point-particle in the absence of radiative corrections has Fi = 1 and F2 = 0, and 
hence Dirac's value for the magnetic moment. The anomalous magnetic moment in Eq. (j59|) 
corresponds to ^2(0) = a/27r. 

One infers from Eq. (|62p that an anomalous contribution to the magnetic moment 
can be associated with an additional interaction term: 



d!^x\q'il){x)a^y'il:{x)F^y{x) , 



(63) 



where F^i,{x) is the gauge-boson field strength tensor. This term is invariant under local 
U{1) gauge transformations but is not generated by minimal substitution in the action for 
a free Dirac field. 

Consider the effect of the global chiral transformation ■0(x) — )• exp(i075)'0(x). The 
term in Eq. (|60p is invariant. However, the interaction of Eq. (|63p is not. These observations 
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facilitate the understanding of a general result: F2 = for a massless fermion in a quantum 
field theory with chiral symmetry realized in the Wigner mode; i.e., when the symmetry 
is not dynamically broken. A firmer conclusion can be drawn. For m = it follows from 
Eq. (|62p that Eq. (|60p does not mix with the helicity-flipping interaction of Eq. (j63|) and 
hence a massless fermion does i iot p ossess a measurable magnetic moment. 



A reconsideration of Ref. 121[ reveals no manifest conflict with these facts. The per- 
turbative expression for F2 (0) contains a multiplicative numerator factor of m and the usual 
analysis of the denominator involves steps that are only valid for m 7^ 0. Fundamentally, 
there is no conundrum because QED is not an asymptotically free theory and hence does 
not have a well-defined nonperturbative chiral limit. 

On the other hand, in QCD the chiral limit is rigorously defined nonperturbatively 



84( 1 ; and the analogue of Schwinger's one- loop calculation can be carried out to find an 
anomalous c/jromo-magnetic moment for the quark. There are two diagrams in this case: 
one similar in fo rm t o that in QED; and another owing to the gluon self-interaction. One 



reads from Ref. 124i | that the perturbative result vanishes in the chiral limit. However, 
Fig.[T] demonstrates that chiral symmetry is dynamically broken in QCD and one must 
therefore ask whether this affects the chromomagnetic moment. 

Of course, it does, and it is now known that this is signalled by the a ppearance of A^s 



in E q. ()56p . Contemporary simulations of lattice-regularized QCD 125| and DSE studies 
|12d | agree that 



Mp,p;0)^-r^B{p^,C) (64) 



and also on the form of Xi{p,p; 0), which is functionally similar to A(p'^, Q. However, owing 
to non-orthogonality of the tensors accompanying Ai and A2, it is difficult to obtain a lattice 
signal for A2- We therefore consider the DSE prediction for A2 in Ref. [126 ] more reliable. 

As pointed out above, perturbative massless-QCD conserves helicity so the quark- 
gluon vertex cannot perturbatively have a term with the helicity-flipping characteristics of 
A3. Equation ()64p is thus remarkable, showing that the dressed-quark-gluon vertex contains 
at least one chirally-asymmetric component whose origin and size owe solely to DCSB; and 
Sec. I VI VI- II illustrates that A3 has a mate rial impact on the hadron spectrum. 

This reasoning is extended in Ref. [123l |: massless fermions in gauge fleld theories 
cannot possess an anomalous chromo/electro-magnetic moment because the term that de- 
scribes it couples left- and right-handed fermions; however, if chiral symmetry is strongly 
broken dynamically, then the fermions should also posses large anomalous magnetic mo- 
ments. Such an effect is expressed in the dressed-quark-gluon vertex via a term 

r^""'(p/,Pi;^) = ^f.ukuT5{pf,p^,k) . (65) 

That QCD generates a strongly momentum-dependent chromomagnetic form facto r 
in the quark-gluon vertex, t^, with a large DCSB-component, is confirmed in Ref. 1251]. 



Only a particular kinematic arrangement was readily accessible in that lattice simulation 
but this is enough to learn that, at the current-quark mass considered: t^ is roughly two 
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orders-of-magnitude larger than the perturbative form; and 
V/>0: \T5{p,-p;2p)\>\X3ip,p;0)\. 



(66) 



The mag nitude of the lattice result is consistent with instanton-liquid model estimates 



127l.ll28l|. 

This large chrom omagneti c moment is likely to have a broad impact on the properties 



of light-quark systems jl28l . Il29l |. In particular, as will be illustrated in Sec. lVI VI-31 it can 



explain the longstanding puzzle of the mass splitting between the ai- and /9-mesons in the 
hadron spectrum [130]. Here a different novel effect will be elucidated; viz., the manner in 
which the quark's chromomagnetic moment generates a quark anomalous e/eciromagnetic 
moment. This demo nstr ation is only possible now that the method of Ref. 83] is available. 
It was accomplished 1231 ] using the same simplification of the effective interaction in Ref. 84 1 
that produced Figs. [71 

In order to understand the vertex Ansatz used in Ref. 123l |. it is necessary to return 
to perturbation theory. One can determine from Ref. 1241 ] that at leading-order in the 
coupling, as, the three-gluon vertex does not contribute to the QCD analogue of Eq. ([59]) 
and the Abelian-like diagram produces the finite and negative correction (— Q;s/[127r]). The 
complete cancellation of ultraviolet divergences occurs only because of the dynamical gen- 
eration of another term in the transverse part of the quark-gluon vertex; namely, 



^t^'^'iPf^Pi) = [^^^ ■ k + i7>^p4A:p]T4(p/,Pi) , 



(67) 



with Tfj,,y = 6f^u - k^ky/k^, a^^ := T^iyU,,. 

Cognisant of this, one may use a simple Ansatz to express the dynamical generation 
of an anomalous chromomagnetic moment via the dressed-quark gluon vertex; viz., 

r.iPf^Pr) = r;^C(p^,p.) + r--(p^,p,), 

with T^{pf,Pi) = (— 7/4)Ab(ppP^), as discussed above, and 



(69) 



u{pf,Pi) =^{z) 



1-2?? 



j^AB{pj,pi) - AA{Pf,pf, 



{p1 + p} 



(70) 



2Ml)/A%, hr = 1 GeV, 



M {s)} is the Euclidean constituent- 



The damping factor F{z) = (1 — exp(— z))/z, z - 
simplifies numerical analysis; and Me = {s\s > 0, s 
quark mass. 

A confined quark does not possess a mass-shell (Sec. lIIip . Hence, one cannot unam- 
biguously assign a single value to its anomalous magnetic moment. One can nonetheless 
compute a magnetic moment distribution. At each value of p'^, spinors can be defined to 
satisfy the free-particle Euclidean Dirac equation with mass m — )• M{p'^) =: <^, so that 



u{pf;^)T^,{pf,Pi;k)u{pi;(;) 



u{pf)[Fi{k^)j^ + — a^ukr.F2{k^)]u{pi) 



(71) 
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FIG. 8: Left panel - /^ (GeV^i) in Eq. ^ cf. {~7/4:)Ab{p'^,p^), both computed using Eq. dH]) and 
the same simphfication of the interaction in Ref. [84] . Right panel - Anomalous chromo- and electro- 
magnetic moment distributions for a dressed-quark, computed using Eq. (j77p . The d ashed -curve in 
both panels is the rainbow-ladder (RL) truncation result. (Figure adapted from Ref. 123|.) 



and then, from Eqs. (|68|) - ([70 
-2^ 7^6% 



H?) 



2?25^^ + 2^6% 



a 



(72) 



where a''^ 



A ^^ ""A 

'' — AAi'i, ?), etc. The numerator's simplicity owes to a premeditated 



cancellation between T4 and T5, which replicates the one at leading-order in perturbation 
theory. Where a compar ison of te rms is possible, this vertex Ansatz is semi-quantitatively 
in agreement with Refs. 1251 . Il26l |. Howe ver, the presence and understanding of the role of 
r^'^™'* is a novel contribution by Ref. [l23l |. NB. It is apparent from Eq. ([72]) that k'^'^ oc m^ 
in the absence of DCSB, so that K^^"^/[2m] — t- in the chiral limit. 



The BSE for the quark-photon vertex can be written following the method of Ref. 82l |. 
Since the method guarantees preservation of the Ward-Takahashi identities, the general 
form of the solution is 






r'^''{Pf,Pi) + rJ;{Pf,P^ 



^„ iPf,Pi) = 7/. -^1 + cr^lukuF2 + T^pcr^ 



pv^v < 



A; F3 + {(l-i ■ k + ij'^a^piukp]Fi 



-ii'^F^ + r^l-kl-kFe 



n- 






n'^l^p- 



Ip '^Vp'' 

^■ukpFs 



(73) 



(74) 



where {Fi\i = 1, . . . , 8} are scalar functions of Lorentz-invariants constructed from pj, pi, 
k. The Ward-Takahashi identity is plainly satisfied; viz.. 



kpi^p{Pf,Pi) = S ^{pf) - S ^{pi) . 



(75) 



Figure [8] depicts the results obtained for the quark's anomalous electromagnetic mo- 
ment form factor 



Mp) ■■-- 



lim -77 

Pf->P 12 K^ 



ii a pykpTl{pf,p) 



h + \p'h ■ 



(76) 
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The result is evidently sizable. It is worth reiterating that /^ is completely nonperturbative: 
in the chiral limit, at any finite order in perturbation theory, /^ = 0. For contrast the 
figure also displays the result computed in the rainbow-ladder truncation of QCD's DSEs. 



As the leading-order in a systematic but stepwise symmetry-preserving scheme [57|, this 
truncation only partially expresses DCSB: it is exhibited by the dressed-quark propagator 
but not present in the quark-gluon vertex. In this case /^ is nonzero but small. These are 
artefacts of the truncation that cannot be remedied at any finite order of the procedure in 



Ref. [53] or a kindred scheme. 

Employing Eq. (|7ip . in connection with the dressed-quark-photon vertex, one can 
write an expression for the quark's anomalous electromagnetic moment distribution 

.M = Mi±f£i±^, pr> 

where: Ak(<?) = 2?^ 5^ — 2?(5^ — "^Fs — ?^-FV; and the Fi are evaluated at ph = pf = 
MijPrY ='■ ?^, k"^ = 0. Plainly, k{<;) = in the chiral limit when chiral symmetry is not 
dynamically broken. Moreover, as a consequence of asymptotic freedom, k(<j) — )• rapidly 
with increasing momentum. The computed distribution is depicted in Fig. [HI It yields 
Euclidean mass-shell values: 

<ii = 0.44 GeV, K-ff = -0.22 , k-- = 0.45 
cf . MIl = 0.35 GeV, Kfg" = , Kfff = 0.048. 

It is thus apparent that DCSB produces a dressed light-quark with a momentum- 
dependent anomalous chromomagnetic moment, which is large at infrared momenta. Signif- 
icant amongst the consequences is the generation of an anomalous electromagnetic moment 
for the dressed light-quark with commensurate size but opposite sign. (NB. This result 
was anticipated in Ref. [I3ll ]. which argued that DCSB usually triggers the generation of a 
measurable anomalous magnetic moment for light-quarks.) The infrared dimension of both 
moments is determined by the Euclidean constituent-quark mass. This is two orders-of- 
magnitude greater than the physical light-quark current-mass, which sets the scale of the 
perturbative result for both these quantities. 

There are two more notable features; namely, the rainbow-ladder truncation, and low- 
order stepwise improvements thereof, underestimate these effects by an order of magnitude; 
and both the t^ and T5 terms in the dressed-quark-gluon vertex are indispensable for a 
realistic description of hadron phenomena. Whilst a simple interaction was used to illustrate 
these outcomes, they are robust. 

These results are stimulating a reanalysis of hadron elastic and transition electromag- 
netic form factors, and the hadron spectrum, results of which will be described below. Fur- 
thermore, given the magnitude of the muon "(/^ — 2 an oma ly" and its assumed importance 
as an harbinger of physics beyond the Standard Model 1321 ] . it might also be worthwhile to 



make a quantitative estimate of the contribution to g^ — 2 from the quark's DCSB-induced 
anomalous moments following, e.g., the computational pattern for the hadro nic l ight-by- 
light scattering component of the photon polarization tensor indicated in Ref. 1331 ] . 



982 SELECTED HIGHLIGHTS FROM THE . . . VOL. 49 



VI-3. ai-p mass splitting 

The analysis in Ref. [84] enables one to construct a symmetry-preserving kernel for the 
BSE given any form for F^. Owing to the importance of symmetries in forming the spectrum 
of a quantum field theory, this is a pivotal advance. One may now use all information 
available, from any reliable source, to construct the best possible vertex Ansatz. The 
last section illustrated that this enables one to incorporate crucial nonperturbative effects, 
which any finite sum of contributions is incapable of capturing, and thereby prove that 
DCSB generates material, momentum-dependent anomalous chromo- and electro-magnetic 
moments for dressed light-quarks. 

The vertex described in Sec. lVI VI-2l contains a great deal of information about DCSB. 
It is the best motivated Ansatz to date and may be used in the calculation of the masses 
of ground-state spin-zero and -one light-quark mesons in order to illuminate the impact of 
DCSB on the hadron spectrum. This analysis expands significantly on the discussion of 
scalar and pseudoscalar mesons in Sec. lVI Vl-li 

A prediction for the spectrum follows once the gap equation's kernel is specified and 
the Ward-Identity solved for the Bethe-Salpeter kernel. A, in Fig. [HI In the axial-vector and 
pseudoscalar channels, the latter is accomplished with A5^(^) = A^S x + Ag^T \, with the 

explicit form of A^9 given in Ref. [84], AgS^ constructed following the procedure described 
therein, and 

Kernels for othe r cha nnels are readily constructed. 

Reference 130l | computes ground-state masses using the method detailed in Ref. 134 |. 



which ensures one need only solve the gap and Bethe-Salpeter equations at spacelike mo- 
menta. This simplifies the numerical problem. To explain, the inhomogeneous BSE is solved 
for the complete Bethe-Salpeter amplitude in a particular channel on a domain of spacelike 
total-momenta, P^ > 0. Any bound-state in that channel appears as a pole in the solution 
at P^ = — m^ggQjj. Denoting the leading Chebyshev moment of the amplitude's dominant 
Dirac structure by T{k; P), where k is the relative momentum, then 1/T[k = 0; P^) exhibits 
a zero at (— w-meson)- The location of that zero is determined via extrapolation of a Fade 
approximant to the spacelike-behavior of l/T{k = 0;P^). This is illustrated for the p- and 
oi-channels in Fig.[9l 

A full set of results is listed in Table HI One first notes the level of agreement between 
Cols. 3 and 4. This illustrates the efficacy of the extrapolation method used to compute 
masses: no diff'erence is greater than 1%. 

The row associated with nia is worth considering next. First compare Cols. 1-3. 
It is an algebraic result that in the RL-truncation of QCD's Dyson-Schwinger equations 
(DSEs), rria ~ 2M, where M is a constituent-like quark mass [138,] . On the other hand, 
incorporating the quark mass function into the Bethe-Salpeter kernel via FJ^^ generates a 
strong spin-orbit interaction, which significantly boosts rrif^ (see Sec. lVI VI-ip . This feature 
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FIG. 9: Illustration of the procedure used to determine meson masses. Solid curve - oi-meson, 
nonperturbative kernel; dot-dash-dash - ai, kernel derived from Eq. ()56p only (Ball-Chiu, BC); and 
dash - oi, leading-order kernel (rainbow-ladder, RL). Dot-dash curve - p-meson, nonperturbative 
kernel; Dot-dash-dot - p, BC-kernel; and dotted - p, RL-kernel. Points - values of l/r(fc = 0; P'^) in 
the given channel computed with the kernel described. Fade approximants are constructe d in e ach 
case; and the location of the zero is identified with {^m^^^^^J. (Figure adapted from Ref. |130| .) 



is evidently unaffected by the inclusion of Tf^"^; i.e., those terms associated with a dressed- 
quark anomalous chromomagnetic moment. Since terms associated with pion fin al-st ate 
interactions are deliberately omitted in the nonperturbative kernel derived in Ref. 1301 ] . it 
is noteworthy that nia in Col. 1 matches estimates for the mass of the dressed- quark-cor e 
component of the o"- meson obtained using unitarised chiral perturbation theory 118l . lll9l |. 
NB. In addition to provi ding a width, such final-state interactions necessarily reduce the 



real part of the mass ^, Il37 1 . 



Now compare the entries in Rows 2, 4-6. The p- and ai -mesons have been known 
members of the spectrum for more than thirty years and are typically judged to be parity- 
partners; i.e., they would be degenerate if chiral symmetry were manifest in QCD. Plainly, 
they are not, being split by more than 400 MeV (i.e., > mp/2). It is suspected that this large 
splitting owes to DCSB. Hitherto, however, no symmetry-preserving treatment of bound- 
states could explain the splitting. This is illustrated by Cols. 3, 4 in the Table, which show 
that whilst a good estimate of mp is readily obtained at leading-order in the systematic 
DSE truncation scheme of Ref. [57|, the axia l -vec tor masses are much underestimated. The 
defect persists at next-to-l eadin g-order liol . Ill2l ]. 

The analysis in Ref. 130l | points to a remedy for this longstanding failure. Using the 
Poincare-covariant, symmetry preserving formulation of the meson bound-state problem 
enabled by Ref. ^83], with nonperturbative kernels for the gap and Bethe-Salpeter equa- 
tions, which incorporate effects of DCSB that are impossible to capture in any step-by-step 
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Rcf. fl30l 


Expt. 


RL-Pade 


RL-direct 


m^ 


0.137 


0.138 


0.138 


0.137 


rup 


0.790 ±0.003 


0.777 


0.754 


0.758 


m„ 


1.08 


0.4- 1.2 


0.645 


0.645 


■niai 


1.27 ±0.02 


1.24 ±0.04 


0.938 


0.927 


mbi 


1.39 ±0.01 


1.21 ±0.02 


0.904 


0.912 


mai - rUp 


0.48 ±0.02 


0.46 ±0.04 


0.18 


0.17 



TABLE I: Computed masses. C ol. 1: Spectrum obtained with the fuU nonperturbative B ethe- 
Salpeter kernels described in Ref. [l30j, which express effects of DCSB. The method of Ref. 134 1 
was used. If noticeable, a jackknife error estimate is reported. Col. 2 - Experimental values; 
computed, except TOo-, from isospin mass-squared averages 1135]. Col. 3 - Masses determined from 
the inhomogeneous BSE at leading-order in the DSE truncation scheme of Re f. 15 TI (with this simple 
kernel, the jackknife error is too small to report); and Col. 4 - results in Ref. [136|, obtained directly 
from the homogeneous BSE at the same order of truncation. 



procedure for improving upon the rainbow-ladder truncation, i t pro vides realistic estimates 
of axial-vector meson masses. In obtaining these results, Ref. 13Cll | showed that the vertex 
Ansatz used most widely in studies of DCSB, r^*-^, is inadequate as a tool in hadron physics. 



Used alone, it increases both nip and ruai but yields rua 



■nir 



0.21 GeV, qualitatively 



unchanged from the rainbow-ladder result (see Fig.[9|). A good description of axial- vector 
mesons is only achieved by including interactions derived from F^*^™; i.e., connected with 



the dressed-quark anomalous chromomagnetic moment [123l |. Moreover, used alone, neither 
term in Y^^ can produce a satisfactory result. The full vertex Ansatz and the associated 
gap and Bethe-Salpeter kernels described in Sec. lVIVI-21 are the minimum required. 

Row 5 contains additional information. The leading-covariant in the 5i-meson channel 
is 75 fc^- The appearance of /c^ suggests immediately that dressed-quark orbital angular 
momentum will play a significant role in this meson's structure, even more so than in 
the oi-channel for which the dominant covariant is 757^. (NB. In a simple quark-model, 
constituent spins are parallel within the oi but antiparallel within the &i. Constituents 
of the 61 may therefore become closer, so that spin-orbit repulsion can exert a greater 
influence.) This expectation is borne out by the following: with the full kernel, ttt-^j increases 
significantly more than rua^ ; and rrib-^ is far more sensitive to the interaction's momentum- 
space range parameter than any other state, decreasing rapidly as the interaction's spatial- 
variation is increasingly suppressed. 

The results reviewed in this section rest on an Ansatz for the quark-gluon vertex 
and whilst the best available information was used in its construction, improvement is 
nonetheless possible. That will involve elucidating the role of Dirac covariants in the quark- 
gluon vertex which have not yet been considered and of resonant contributions; viz., meson 
loop effects that give widths to some of the states considered. In cases for which empirical 
width-to-mass ratios are < 25%, one might ju dge t hat such contributions can reliably be 
obtained via bound-state perturbation theory [1391 ] . Contemporary studies indicate that 
these effects reduce bound-state masses but the reduction can uniformly be compensated 
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by a modest inflation of the interaction's mass-scale [l38|, |l40(] , so that the masses in Table |T] 
are semiquantitatively unchanged. The case of the cr-meson is more complicated. However, 
the prediction of a large mass for this meson's dre s sed- quark core leaves sufficient room for 



a strong reduction by resonant contributions 118l . Ill9l | . 



This section reviewed a continuum framework for computing and explaining the meson 
spectrum, which combines a veracious description of pion properties with estimates for 
masses of light-quark mesons heavier than nip- (A contemporary lattice-QCD perspective 
on this problem may be drawn from Ref. [141i] .) The method therefore offers the promise 
of a first reliable Poincare- invariant, symmetry-preserving computation of the spectrum of 
light-quark hybrids and exotics; i.e., those putative states which are impossible to construct 
in a quantum mechanics based upon constituent-quark degrees-of-freedom. So long as the 
promise is promptly fulfilled, the approach will provide predictions to guide the forthcoming 
generation of facilities. 



VII. PION ELECTROMAGNETIC FORM FACTORS 

In charting the long-range interaction between light-quarks via the feedback between 
experiment and theory, hadron elastic and transition form factors can provide unique in- 
formation, beyond that obtained through studies of the hadron spectrum. 

VII-1. Charged pion 

This is demonstrated very clearly by an analysis of the electromagnetic pion form 
factor, F^™((5^), because the pion has a unique place in the Standard Model. It is a bound- 
state of a dressed-quark and -antiquark, and also that almost-massless collective excitation 
which is the Goldstone mode arising from the dynamical breaking of chiral symmetry. 
This dichotomy can only be understood through the symmetry-preserving analysis of two- 
body bound-states [M]. Furthermore, the possibility that this dichotomous nature could 
have wide-ranging effects on pion properties has made the empirical investigation of these 
properties highly desirable, despite the difficulty in preparing a system that can act as a 
pioii targ et and the concomitant complexities in the interpretation of the experiments; e.g., 
14i - [l45| . 



The merit of using F^^{Q'^) to elucidate the potential of an interplay between exper- 
iment and nonperturbative theory as a means of constraining the long-range behaviour of 
QCD's /3-function is amplified by the existence of a prediction that Q'^Ft^{Q'^) ~ constant 
for Q^ ^ mj in a theory whose interaction is mediated by massless vector-bosons. To be 
explicit [uillii ]: 



Q^F^iQ^) "^ -'"^ 16vr/>(Q2), (80) 

which takes the value 0.13 GeV^ at Q^ = lOGeV^ if one uses the one-loop result 
q(Q^ = lOGeV ) ^ 0.3 . The verification of this prediction is a strong motivation for 



modern experiment 142Hl44i |. which can also be viewed as an attempt to constrain and 
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TABLE II: Results obtained with (in GeV) tog = 0.132, Aj^ = 0.24, A„ 
quantities are listed in GeV. 



— 0.905. Dimensioned 



m 


E^ 


F^ 


Ep 


M 


1/3 


TO^ 


TOp 


u 


fp 



0.007 


3.568 
3.639 


0.459 
0.481 


1.520 
1.531 


0.358 
0.368 


0.241 
0.243 



0.140 


0.919 
0.928 


0.100 
0.101 


0.130 
0.129 



map experimentally the pointwise behaviour of the exchange interaction that binds the 
pion. 

Section [VV-ll details some extraordinary consequences of DCSB, amongst them the 
Goldberger-Treiman relations of Eqs. JHM - (fT7|) . Of these, Eqs. (fTSl) and (fT6]) entail that 
the pion possesses components of pseudovector origin which alter the asymptotic form of 
-F^™(Q'^) by a multiplicative factor of Q^ cf. the result obtained in their absence [66(. 

QCD-based DSE calculations of F^^(0^ exist 66|, |7j], the most systematic of which 
[TJ] predicted the measured form factor 142] ■ Germane to this discourse, however, is 
an elucidation of the sensitivity of -F^™(Q^) to the pointwise behaviour of the interaction 
between quarks. We therefore recapitulate on Refs. 68|,l70|, which explored how predictions 
for pion properties change if quarks interact not via massless vector-boson exchange but 
instead through a contact interaction; viz., 



1 



9 Df,„{p-q) = 5f,^^ , 



m 



(81) 



G 



where mo is a gluon mas s-scale (such a scale is generated dynamically in QCD, with a 
value ~ 0.5 GeV [481 . |l49t]), and proceeded by embedding this interaction in a rainbow- 
ladder truncation of the DSEs. 

In this case, using a confining regularisation scheme |l50l]. the gap equation, which 
determines this interaction's momentum-independent dressed-quark mass, can be written 



M = m + 



M 



3^2 



m 



■C{M^;n„T^ 



G 



where m is the current-quark mass and 

C(M2;rir,ruv) = M2[r(-1,M2T2J -r(-l,MV2)] 



(82) 



(83) 



with T{a,y) being the incomplete gamma-function. Results are presented in Tablellll 

With a symmetry-preserving regularisation of the interaction in Eq. (j8ip , the Bethe- 
Salpeter amplitude cannot depend on relative momentum. Hence Eq. (|13p becomes 



r.(P) = 75 



iE^{P) + —l-PF^iP) 



(84) 
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and the explicit form of the model's ladder BSE is 

1 



FAP) 



37r2 



rri] 



G 



K-EE ^EF 
K-FE K,FF 



E^{P) 

PAP) 



where, with m = = P^, anticipating the Goldstone character of the pion, 



2}Cfe = Ci{M^;4,tI), JCff 



-2ICfe , 



(85) 



(86) 



and Ci{z) = —zC'{z), where we have suppressed the dependence on Tji-^uv- The solution of 
Eq. (j85p gives the pion's chiral-limit Bethe-Salpeter amplitude, which, for the computation 
of observables, should be normalised canonically; viz.. 



Pa=N,tr 



(27r) 
Hence, in the chiral limit, 

Nr 1 



^^A-P)^Siq + P) KiP) Siq) . 



/« 



1 



Ci{M';t-^,t^,)E4E^ - 2F^], 



47r2M2 
and the pion's leptonic decay constant is 

Nr 1 



(87) 



(88) 



J w 



47r2M 



CiiM'-T^^,Ti^)[E^-2F^] 



If one has preserved Eq. ([8]), then, for rn, = in the neighbourhood of P^ = 0, the 
solution of the axial-vector BSE has the form: 



P 



and the following subset of Eqs. ()14p - (J17p will hold: 



f^E^ = M, 2-^+Fr = 1. 

-C/TT 



(90) 



(91) 



Hence F-k{P) is necessarily nonzero in a vector exchange theory, irrespective of the pointwise 
behaviour of the interaction. It has a measurable impact on the value of f-,^. 

Based upon these results, one can proceed to compute the electromagnetic pion form 
factor in the generalised impulse approximat ion 68|,|70|]; i.e., at leading-order in a symmetry- 
preserving DSE truncation scheme |66l . l74l . Il5ll |. Namely, for an incoming pion with mo- 
mentum pi = K — Q/2, which absorbs a photon with space-like momentum Q, so that the 
outgoing pion has momentum p2 = K -\- Q/2, 



2K^F^^{Q^) = 2N, 



dH 



-tro 



ir^{-p2)S{t+p2)i^^,PT{Q'')S{t+pi) iV^ipi) Sit) 



, (92) 
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FIG. 10: Dressing function for the transverse piece of the quark-photon vertex, computed using the 
parameter values described in connection with Table [Hi The pole associated with the ground-state 
vector meson is clear. (Figure adapted from Ref. |70|.) 



where Pt(Q^)) depicted in Fig.[TOl describes the full extent of dressing on the quark-photon 
vertex produced by a contact interaction in the rainbow-ladder truncation [TC?]. The form 
factor is expressible as follows: 

^2\ 



FZ^iQ') 



FZee{Q^) + FZef{Q^) + FZff{Q') 



(93) 
(94) 
(95) 



where F2"'^{Q'^ 



^,^^ 



is that part of the form factor produced by the undressed quark-photon 
vertex and the functions T'^ have simple algebraic forms in this calculation [TO]. 

In the left panel of Fig. [11] we present F^"^{Q'^) and the three separate contributions 
defined in Eq. (I95p . It is evident that, in magnitude, F!^"j^p contributes roughly one-third of 
the pion's unit charge. This could have been anticipated from Eq. (j88p . More dramatically, 
perhaps: the interaction in Eq. (fSTj) generates 



F^iQ^ ^ oo) = constant. 



(96) 



Both results originate in the nonzero value of F^^ (P) , which is a straightforward consequence 
of the symmetry-preserving treatment of a vector exchange theory [64] . Equation (j96p 
should not come as a surprise: with a symmetry-preserving regularisation of the interaction 
in Eq. (jSTj) . the pion's Bethe-Salpeter amplitude cannot depend on the constituent's relative 
momentum. This is characteristic of a pointlike particle, which must have a hard form 
factor. The right panel of the figure illustrates that the necessary inclusion of P't{Q^) is 
critical in the timelike region and has a measurable quantitative impact for a significant 
range of spacelike momenta. It does not, however, affect the truly ultraviolet behaviour. 
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FIG. 11: Left panel. F^'^{Q^) and the separate contributions introduced in Eq. dM])- F^'^{Q^ = 
0) = 1, without fine-tuning, because a symmetry-preserving regularisation of the interaction in 
Eq. (|8ip was implemented. Right panel. -F°™((5^) computed in rainbow-ladder truncation from the 
interaction in Eq. (j8ip : solid curve - fully consistent, i.e., with a dressed-quark-photon vertex so 
that the p-pole appears; and dashed curve - computed using a bare quark-photon vertex, namely, 
F^"j{Q'^). Dotted curve - fit to the result in Ref. [745, which also included a consistently-dressed 
quark-photon vertex and serves to illustrate the trend of contemporary data. (Figures adapted from 
Refs. [silel.) 



In Fig. 1121 we compare the form factor computed from Eg. (I81[) with contemporary 
experimental data 142Hl44i | and a QCD-based DSE prediction (74] . Both the QCD-based 
result and that obtained from the momentum-independent interaction yield the same values 
for the pion's static properties |68l47Cll| . However, for Q^ > the form factor computed 
using ~ l/A;^-vector-boson exchange is immediately distinguishable empirically from that 
produced by a momentum-independent interaction. Indeed, the figure shows that for -F^™, 
existing experiments can already distinguish between different possibilities for the quark- 
quark interaction. 

Combining Figs.llllandll2lit becomes apparent that F^"ee ^^ only a good approxima- 
tion to the net pion form factor for Q^ < M^. F^'^^ and F^^p evolve with equal rapidity - 
there is no reason for this to be otherwise, as they are determined by the same mass-scales 
- but a nonzero constant comes quickly to dominate over a form factor that falls swiftly to 
zero. 

It is plain now that when a momentum-independent vector-exchange interaction is 
regularised in a symmetry-preserving manner, the results are directly comparable with ex- 
periment, computations based on well-defined and systematically-improvable truncations 
of QCD's DSEs TJ], and perturbative QCD. In this context it will be apparent that 



a contact interaction, whilst capable of describing pion static properties well, TableHU 
generates a form factor whose evolution with Q^ deviates markedly from experiment for 
Q^ > 0.17 GeV^ ^ M^ and produces asymptotic power-law behaviour, Eq. (f96]l 



m serious 



conflict with QCD |146l4l48l |. 

In that connection Fig. [T] and Eqs. (114p - (I17p are relevant again. In the electromag- 
netic elastic scattering process, the momentum transfer, Q, is primarily shared equally be- 
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Q^ [GeV^] 



FIG. 12: Solid curve : Q^F^j{Q^) obtained with Eq. dH]). Dashed curve : DSE prediction [74.], which 
employed a momentum-dependent renormahsation-group-improved gluon exchange interaction. For 
Q^ > 0.17 GeV^ « M^, marked by the vertical dotted line, the contact int eraction r esult for -F°™((5^) 
differs from that in Ref. 
fromRef. S.) 



[ll] by more than 20%. The data are from Refs. [l42l4l44J . (Figure adapted 



tween the pion's constituents because the bound-state Bethe-Salpeter amplitude is peaked 
at zero relative momentum. Thus, one can consider k ~ Q/2. The Goldberger-Trieman-like 
relations express a mapping between the relative momentum of the pion's constituents and 
the one-body momentum of dressed-quark; and the momentum dependence of the dressed- 
quark mass function is well-described by perturbation theory when k'^ > 2GeV^. Hence, 
one should expect a perturbative-QCD analysis of the pion form factor to be valid for 



k' = Q'/4>2Ge\';le., 



(97) 



This explains the result in Ref. [66|]. A similar argument for baryons suggests that the 
nucleon form factors should be perturbative for Q^ > 18GeV^. 

It is worth reiterating that the contact interaction produces a momentum-independent 
dressed-quark mass function, in contrast to QCD-based DSE studies [^jlldl and lattice-QCD 
141 ] . This is fundamentally the origin of the marked discrepancy between the form factor it 
produces and extant experiment. Hence Refs. [6814701] highlight that form factor observables, 
measured at an upgraded Jefferson laboratory, e.g., are capable of mapping the running 
of the dressed-quark mass function. Efforts are underway to establish the signals of the 
running mass in baryon elastic and transition form factors. 
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VII-2. Neutral pion 

The process 7*7 — t- vr'^ is also of great interest because in order to explain the asso- 
ciated transition form factor within the Standard Model on the full domain of momentum 
transfer, one must combine, using a single internally-consistent framework, an explanation 
of the essentially nonperturbative Abelian anomaly with the features of perturbative QCD. 
The case for attempting t his r eceived a significant boost with the publication of data from 



the BaBar Collaboration 154l | because, while the y agree w ith earlier experiments on their 



common domain of squared- momentum-transfer 155l . Il56l | , the BaBar data are unexpect- 
edly far above the prediction of perturbative QCD at larger values of Q^. 



This so-called "Babar anomaly" was considered in Ref. 69(], wherein it is argued that 



in fully-self-consistent treatments of pion: static properties; and elastic and transition form 
factors, the asymptotic limit of the product Q'^G^*^t^o{Q'^), which is determined a priori 
by the interaction employed, is not exceeded at any finite value of spacelike momentum 
transfer: the product is a monotonically-increasing concave function. A consistent approach 
is one in which: a given quark-quark scattering kernel is specified and solved in a well- 
defined, symmetry-preserving truncation scheme; the interaction's parameter(s) are fixed 
by requiring a uniformly good description of the pion's static properties; and relationships 
between computed quantities are faithfully maintained. 

Within such an approach it is nevertheless possible for Q'^F^'^{Q'^) to exceed its 
perturbative- QCD asymptotic limit because the leading-order matrix-element involves two 
Bethe-Salpeter amplitudes. This permits an interference between dynamically-generated 
infrared mass-scales in the computation. Moreover, for F"^^(Q^) the perturbative QCD 
limit is more than an order-of-magnitude smaller than rni. Owing to the proximity of the 
p-meson pole to Q^ = 0, the latter mass-scale must provide a fair first-estimate for the 
small-Q^ evolution of F^™((5^). A monopole based on this mass-scale exceeds the pQCD 
limit VQ^ > 0. For the transiti on fo rm factor, however, the opposite is true because rni is 
less-than the pQCD limit; viz. 1481 ] . 



lim Q'2GiQ') = 87r'f^. (98) 

The vector current-current contact-interaction canvassed in this Section may be de- 
scribed as a vector-boson exchange theory with vector-field propagator (1/A:^)'', k = 0. It 
was shown [69|] that the consistent treatment of such an interaction produces a 7*7 — t- vr'' 
transition form factor that disagrees with all available data. On the other hand, precisely 
the same treatment of an interaction which preserves the one-loop renormalisation group 
behaviour of QCD, produces a for m fac tor in good agreement with all but the large-Q^ 



data from the BaBar Collaboration 154i |. These points are illustrated in Fig. 1131 

Studies exist which interpret the BaBar data as an indication that the pion's dis- 
tribution amplitude, (pnix), deviates dramatically from its QCD asympt otic form, indeed. 



that 4>Tr{x) = constant, or is at least flat and nonvanishing at x = 0, 1 [1581 . Il59l |. How- 
ever, it has often been explained |9l. I67H7Q] that such a distribution amplitude characterises 
an essentially-pointlike pion; and, as we have seen, when used in a fully-consistent treat- 
ment, it produces results for pion elastic and transition form factors that are in striking 
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FIG. 13: Q -weighted 7*7 — >■ n transition form factor. Data: red circles, Ref. [15j|; green di- 
amonds, Ref. [155[; and blue squares, Ref. |156| . Solid curve - Q'^G{Q'^) computed using the 
symmetry-preserving, fuhy-self-consistent rainbow-ladder treatment of the contact interaction in 
Eq. (|81l) . which produces a pion distribution ampUtude (j)Tr{x) = constant; and dot-dashed curve - fit 
to the 7*7 — > TT*^ transition form factor computed in a QCD-based rainbow-ladder-truncation DSE 
study [157]. Both curves have been divided by {2'k'^ f^) in order to match the data's normalisation. 
(Figure adapted from Ref. t69|.) 



disagreement with experiment. Reiterating, a bound-state pion with a pointlike component 
will produce the hardest possible form factors; i.e., form factors which become constant at 
large-Q^. 

On the other hand, QCD-based studies produce soft pions, a valence-quark distribu- 
tion amplitude for the pion that vanishes as ~ (1 — x)^ for x ~ 1, and results that agree 
well with the bulk of exist iiig data. We will return to this in the next section. 

The analysis in Ref. [6{|] shows that the large-Q^ BaBar data is inconsistent with 
QCD and also inconsistent with a vector current-current contact interaction. It supports 
a conclusion that the large-Q^ data reported by BaBar is not a true repr esentatio n of the 



7 7 



TT" 



transition form factor, a perspective also developed elsewhere |16CHl63l |. There 



ri'l 



is experimental evidence in support of this view; namely, the 7* — )■ r/7 and 7* 
transition form factors have also been measured by the BaBar Collaboration (164], at Q^ 
112 GeV^, and in these cases the results from CLEO [156] and BaBar are fully consistent 
with perturbative-QCD expectations. 
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VIII. PION AND KAON VALENCE-QUARK DISTRIBUTIONS 

The past forty years have seen a tremendous effort to deduce the parton distribution 
functions (PDFs) of the most accessible hadrons - the proton, neutron and pion. There are 
many reasons for this long sustained and thriving interest [9] but in large part it is moti- 
vated by the suspected process-independence of the usual parton distribution functions and 
hence an ability to unify many hadronic processes through their computation. In connection 
with uncovering the essence of the strong interaction, the behaviour of the valence-quark 
distribution functions at large Bjorken-x is most relevant. Furthermore, an accurate deter- 
mination of the behavior of distribution functions in the valence region is also important 
to high-energy physics. Particle discovery experiments and Standard Model tests with col- 
liders are only possible if the QCD background is completely understood. QCD evolution, 
apparent in the so-called scaling violations by parton distribution functions, entails that 
with increasing center-of-mass energy, the support at large-x in the distributions evolves to 
small-x and thereby contributes materially to the collider background. Signficantly, in the 
infinite momentum frame, Bjorken-x measures the fraction of a hadron's four-momentum 
carried by the struck parton and, e.g., the valence-quark distribution function, qv{x), mea- 
sures the number-density of valence-quarks with momentum-fraction x. NB. The nucleon 
PDFs are now fairly well determined for x < 0.8 but the pion and kaon PDFs remain poorly 
known on the entire domain of x. 

Owing to the dichotomous nature of Goldstone bosons, understanding the valence- 
quark distribution functions in the pion and kaon is of great importance. Moreover, given 
the large value of the ratio of s-to-u current-quark masses, a comparison between the pion 
and kaon structure functions offers the chance to chart effects of explicit chiral sym metry 



breaking on the structure of would-be Goldstone modes. There is also the prediction [153, 
153| that a theory in which the quarks interact via l//c^-vector-boson exchange will produce 



valence-quark distribution functions for which 

gv(a;) oc (l-x)2+^, x >0.85, (99) 

where 7 > is an anomalous dimension that grows with increasing momentum transfer. 
(See Sec.VLB.3 of Ref. [9j for a detailed discussion.) 

Experimental knowledge of the parton structure of the pion and kaon arises primarily 
from pionic or kaonic Drell-Yan processes involving nucleons in heavy nuclei 145l . Il65 - 



1681 ]. Theoretically, given that DCSB plays a crucial role in connection with pseudoscalar 
mesons, one must employ an approach that realistically expresses this phenomenon. The 
DSEs therefore provide a natural framework: studies of the pion and kaon exist and will 



be reviewed here. The first 67| computed pion PDFs, using efficacio us pararnetris ations 



of both the Bethe-Salpeter amplitude and dressed-quark propagators [ l5ll . Il69l . Il70l |. The 



second 7l| employed direct, numerical DSE solutions in the computation of the pion and 
kaon PDFs, ada ptin g t he ap proach employed in successful predictions of electromagnetic 
form factors 59l. I73l. l74l . ll57l |: and also studied the ratio uk{x)/utj{x) in order to elucidate 



aspects of the influence of an hadronic environment. 
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FIG. 14: Pion valence quark distribution function evolved to (5.2 GeV)^. Solid curve - full DSE 
calculation 7l|; dot-dashed curve - se mi-p henomenological DSE-based calculation in Ref. 67|; filled 
circles - experimental data from Ref. 168], at scale (4.05 GeV)^; dashed curve - NLO reanalysis of 
the experimental data 145|; and dot - dot- dashed curve ~ NLO reanalysis of experimental data with 
inclusion of soft-gluon resummation [171| . (Adapted from Ref. [7l|.) 



In rainbow-ladder truncation, one obtains the pion's valence-quark distribution from 



Ut,{x) 






.tr,d [r.(£, -P) Suii) T^ii; x) Su{i) T^{i, P) Sdii - P)] , (100) 



wherein the Bethe-Salpeter amplitudes and dressed-quark propagators are discussed above 
and T^{i;x) is a generalization of the dressed-quark-photon vertex, describing a dressed- 
quark scattering from a zero momentum photon. It satisfies a BSE (here with a rainbow- 
ladder kernel) with the inhomogeneous term i^ ■ n6{i ■ n — xP ■ n). Here n^ is a light-like 
vector satisfying v? = 0. In choosing rainbow-ladder truncation one implements a precise 
parallel to the symmetry-preserving treatment of the pion charge form factor at Q^ = 0, 
wherein the vector current is conserved by use of ladder dynamics at all three vertices and 
rainbow dynamics for all three quark propagators 66|, |73|, |7J, Il5ll | . Equation (llOOp ensures 
automatically that 



(4) := 



dx q'j (a 



1 for / = u, d , 



(101) 



since JdxT^{i;x) gives the Ward-identity vertex and the Bethe-Salpeter amplitudes are 
canonically normalised. 
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Figure [T^ displays the DSE results for the pion's valence u-quark distribution, evolved 
from a resolving scale Qg = (0.6 GeV)^ to Q^ = (5.2 GeV)^ using leading-order evolution, 
and a comparison with vrA^ Drell-Yan data [168!] at a scale Q^ ~ (4.05 GeV)^, inferred via 
a leading-order analysis. The computation's resolving scale, Qo, was fi xed b y matching the 



(x"')'^ moments for n = 1,2,3 to an experimental analysis at (2GeV)^ 1721 ]. 
It is notable that at Qo the DSE results yield 

2(x)5„=0.7, 2(x)^„=0.8. (102) 



(For comparison, the parametrised valence-like pion parton distributions of Ref. [1731 ] yield 
a gluon momentum fraction of {xg)^ =o 5i ~ 0.3.) In each case the remainder of the 
hadron's momentum is carried by gluons, which effect binding of the meson bound state 
and are invisible to the electromagnetic probe. Some fraction of the hadron's momentum 
is carried by gluons at all resolving scales unless the hadron is a point particle [91]. In- 
deed, it is a simple algebraic exercise to demonstrate that the only non-increasing, convex 
function which can produce {x^) = 1 and (x) = ^, is the distribution u{x) = 1, which is 
uniquely connected with a pointlike meson; viz., a meson whose Bethe-Salpeter amplitude 
is momentum- independent. Thus Eqs. (|102p are an essential consequence of momentum 
conservation. 

Whilst the DSE results in Fig.[ll]are both consistent with Eq. ([99]) : i.e., they produce 
algebraically the precise behaviour predicted by perturbative QCD, on the valence-quark 
dom ain it is evident that they disagree markedly with the Drell-Yan data reported in 



Ref. 1681 ]. This tension was long seen as a crucial mystery for a QCD description of the 



lightest and subtlest ha dron [9|. Its re-emergence with Ref. [67|] motivated a NLO reanalysis 
of the Drell-Yan data [1451], the result of which is also displayed in Fig.[T31 At NLO the 
extr acted PDF is softer at high-x but the discrepancy nevertheless remains. To be precise. 



Ref. 1451 ] determined a high-x exponent of /3 ~ 1.5 whereas the exponents produced by the 



DSE studies [67l. l7ll] are, respectively, 2.1 and 2.4 at the common model scale. They do not 
allow much room for a harder PDF at high-x. 

Following the highlighting of this discrepancy in Ref. [9[] , a r esolu tion of the conflict 
between data and well-constrained theory was proposed. In Ref. [171] a long-overlooked 
effect was incorporated; namely, "soft-gluon resummation." With the inclusion of this 
next-to-leading-logarithmic threshold resummation effect in the calculation of the Drell- 
Yan cross section, a considerably softer valence-quark distribution was obtained at high-x. 
This is readily understood. The Drell-Yan cross-section factorises into two pieces: one 
hard and the other soft. The soft piece involves the PDF and the hard piece is calculable 
in perturbation theory. Adding additional interactions to the latter, which are important 
at large-x; viz., soft gluons, provides greater strength in the hard piece on the valence- 
quark domain. Hence a description of the data is obtained with a softer PDF. Indeed, the 
distribution obtained thereby matches precisely the expectations based on perturbative- 
QCD and obtained using DSEs. This is evident in a comparison between the dash-dot and 
dash- dot- dot curves in Fig. (HI 

The ratio uk/ut^ measures the effect of the local hadronic environment. In the kaon, 
the n-quark is bound with a heavier partner than in the pion and this should cause uk{x) 
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FIG. 15: DSE prediction for the ratio of w-quark distributions in the kaon and pion [Ol [DJ. The 
full Bethe-Salpeter amplitude produces the solid curve; a reduced BSE vertex produces the dashed 
curve. The reduced ampHtude retains only the invariants and amphtudes involving pseudoscalar 
and axial vector Dirac matrices, and ignores dependence on the variable q ■ P. These are part of 
the red uctions th at occur in a pointlike treatment of pseudoscalar mesons. The experimental data 
is from [leilie^. (Adapted from Ref. Tl^.) 



to peak at lower-x than Ujix) . Th e DSE prediction [y, f7ll] is shown in Fig. 1151 along with 
available Drell-Yan data 1651 . \l6u!\ . The parameter-free DSE result agrees well with the 
data. We note that 



uxiO) DGLAP:i 



(103) 



viz, the ratio approaches one under evolution to larger r esolv ing scales owing to the increas- 
ingly large population of sea-quarks produced thereby [17j]. On the other hand, the value 
at X = 1 is a fixed-point under evolution: 



vg? > Ql 



ui<'(i) 



UTvil) 



DGLAP:Q2^Q2 Uk{1) 



Qi 



U^{1) 



uk{1) 



Qg 



M7r(l) 



(104) 



Qg 



i.e., it is the same at every value of the resolving scale Q^, and is therefore a persistent 
probe of nonperturbative dynamics [9[. 

With Ref. 7l|] a significant milestone was achieved; viz., unification of the computa- 
tion of distribution functions that arise in analyses of deep inelastic scattering with that 
of numer ous other properties of pseudoscalar mesons, including meson-meson scattering 
175l . Il76l | and the successful prediction of electromagnetic elastic and transition form fac- 
tors. The results confirm the large-x behavior of distribution functions predicted by the 
QCD parton model; provide a good account of the vr-A^ Drell-Yan data for n^(x); and a 
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parameter- free prediction for the ratio uk{x) / Ut^{x) that agrees with extant data, showing 
a strong environment-dependence of the tt-quark distribution. The new Dreh-Yan exper- 
iment running at FNAL is capable of vahdating this comparison, as is the COMPASS II 
experiment at CERN. Such an experiment should be done so that complete understanding 
of QCD's Goldstone modes can be claimed. 



IX. BARYON PROPERTIES 

While a symmetry-preserving description of mesons is essential, it is only part of 
the physics that nonperturbative QCD must describe because Nature also presents us with 
baryons: light-quarks in three-particle composites. An explanation of the spectrum of 
baryons and the nature of interactions between them is basic to understanding the Standard 
Model. The present and planned experimental programmes at JLab, and other facilities 
worldwide, are critical elements in this effort. 

No approach to QCD is comprehensive if it cannot provide a unified explanation of 
both mesons and baryons. We have explained that DCSB is a keystone of the Standard 
Model, which is evident in the momentum-dependence of the dressed-quark mass function 
- Fig.[Tl it is just as important to baryons as it is to mesons. The DSEs furnish the only 
extant framework that can simultaneously connect both meson and baryon observables with 
this basic feature of QCD, having provided, e.g., a direct correlation of meson and baryon 
properties via a single interaction kernel, which prese rves QCD's on e-loop renormalisation 



group behaviour and can systematically be improved [14d . ll77H179l |. 

In quantum field theory a baryon appears as a pole in a six-point quark Green func- 
tion. The residue is proportional to the baryon's Faddeev amplitude, which is obtained from 
a Poincare covariant Faddeev equation that sums all possible exchanges and interactions 
that can take place between three dressed-quarks. A tractable Faddeev equation for baryons 
jl80.] is founded on the observation that an interaction which describes colour-singlet mesons 
also gen erate s nonpointlike quark-quark (diquark) correlations in the colour-3 (antitriplet) 
channel |18ll |. The dominant correlations for ground state octet and decuplet baryons are 
scalar (0^) and axial- vector (l"*") diquar ks b e caus e, for example, the associated mass-scales 



are smaller than the baryons' masses 109l . Il82l | and their parity matches that of these 



baryons. It follows that only they need be retained in approximat ing the q uark-quark 
scattering matrix which appears as part of the Faddeev equation 29|, Il38l . Il77 [. 



We note that diquarks do not appear in the strong interaction spectrum [57l. ll07l . ll08l | 
but the attraction between quarks in this channel justifies a picture of baryons in which 
two quarks within a baryon are always correlated as a colour-3 diquark pseudoparticle, and 
binding is effected by the iterated exchange of roles between the bystander and diquark- 
participant quarks. Here it is in ipor tant to emphasise strongly that QCD supports non- 



pointlike diquark correlations 7d . ll83l |. Hence models that employ pointlike diquark degrees 



of freedom have little connection with QCD. 

Numerous properties of the nucleon have recently been described using the Faddeev 
equation just outlined, e.g.: a survey of the form factors [29|, with additional results expli- 
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cated in Refs. 38l . Il37l . ll84l4l88l | ; in dications of the role of quark orbital angular momentum 



in forming the nucleon's spin [32]; and computation of the nucleon's cr-term jlS]. This 
body of work will be reviewed elsewhere. However, muc h mo re can and should be done, as 



described, for example, in Sec. Ill of Ref. ^189]; and Ref. [138l | is a first step. 

It is expected that feedback between DSE results and extant and forthcoming preci- 
sion data on nucleon elastic and transition form factors can serve as a practical means by 
which to empirically chart the momentum evolution of the dressed-quark mass function, 
and therefrom the infrared behavior of QCD's /3-function. In particular, it should enable 
the unambiguous location of the transition boundary between the constituent- and current- 
quark domains, signalled by the sharp drop with increasing momentum that is apparent in 
Fig.[Tl which can likely be related to an infrared inflexion point in QCD's running coupling, 
whose properties are determined by the /3-function. 



X. EPILOGUE 

Dynamical chiral symmetry breaking (DCSB) is a fact in QCD. It is manifest in 
dressed-propagators and vertices, and, amongst other things, it is responsible for: the 
transformation of the light current-quarks in QCD's Lagrangian into heavy constituent-like 
quark's, in terms of which order was first brought to the hadron spectrum; the unnatu- 
rally small values of the masses of light-quark pseudoscalar mesons; the unnaturally strong 
coupling of pseudoscalar mesons to light-quarks - g.,rqq ~ 4.3; and the unnaturally strong 
coupling of pseudoscalar mesons to the lightest baryons - 5^^^ ss 12.8 ss "ig-wqq- 

Herein we have used a diverse range of phenomena to illustrate the dramatic impact 
that DCSB has upon meson observables, amongst them: the spectrum; form factors; and 
parton distribution functions. Unavoidably, many valuable contributions during th e last 
few years were overlooked, some of which may be pursued through Refs. 190|, ll9l|. We 



have also completely neglected an extensive body of work that focuses on elucidating the 



nat ure of th e primordial state of matter that has been recreated with RHIC: Refs. 25- 



27, 191 



192] serve as a starting point for exploring this area of DSE research. (As noted in 



Ref. [1951 ]. the rainbow truncation of QCD's gap equation excludes what in a point-meson 
field theory would appear as l/A'^c-suppressed meson-loop dressing of the quark propagator, 
and this is why chiral symme try r estoration is a mean field transition in all models within 



the rainbow-truncation class 196l |. Such meson-loop corrections are capable of correcting 



the critical exponents; i.e., of introducing that infrared divergence in the dressed-quark self- 
energy which can force a deviation from mean fie ld cr itical exponents in the chiral symmetry 
restoring transition. This is illustrated in Ref. 197|.) In addition, there are attempts to 



derive an equation of state that may be used in modern astrophy sics i n or der, e.g., to assist 
in the identification of a neutron star with a quark-matter core 193l . Il94 ]. 



A "smoking gun" for DCSB is the behaviour of the dressed-quark mass function. 
The momentum dependence manifest in Fig. [T] is an essentially quantum field theoretical 
effect. Exposing and elucidating its consequences therefore requires a nonperturbative and 
symmetry-preserving approach, where the latter means preserving Poincare covariance, chi- 



VOL. 49 LEI CHANG, CRAIG D. ROBERTS, AND PETER C. TANDY 999 



ral and electromagnetic current-conservation, etc. The Dyson- Schwinger equations (DSEs) 
provide such a framework. Experimental and theoretical studies are underway that will 
identify observable signals of Mijp') and thereby confirm and explain the mechanism re- 
sponsible for the vast bulk of visible mass in the Universe. 

There are many reasons why this is an exciting time in hadron physics. We have 
focused on one. Namely, through the DSEs, we are positioned to unify phenomena as 
apparently diverse as: the hadron spectrum; hadron elastic and transition form factors, 
from small- to large-Q^; and parton distribution functions. The key is an understanding 
of both the fundamental origin of nuclear mass and the far-reaching consequences of the 
mechanism responsible; namely, DCSB. These things might lead us to an explanation of 
confinement, the phenomenon that makes nonperturbative QCD the most interesting piece 
of the Standard Model. 
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